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Abstract 

We present a short review of our basic understanding of the physics of copper 



S 

5;^ ■ oxide superconductors and formulate the hst of "solved" and "unsolved" problems. 

O I The main problem remains theoretical description of the properties of the normal 

state, requiring clarification of the nature of the so called pseudogap state. 

We review simplified models of the pseudogap state, based on the scenario of 
^ , strong electron scattering by (pseudogap) fluctuations of "dielectric" (AFM, CDW) 

0^ I short - range order and the concept of "hot" spots (patches) on the Fermi surface. 

^ • Pseudogap fluctuations are described as appropriate static Gaussian random field 

00 ■ scattering electrons. 

I We derive the system of recurrence equations for the one - particle Green's 

' function and vertex parts, taking into account all Feynman diagrams for electron 

scattering by pseudogap fluctuations. Results of calculations of spectral density, 
j3 ! density of states and optical conductivity are presented, demonstrating both pseu- 

dogap and localization effects. 

■ We analyze the anomalies of superconducting state (both s— and d - wave pair- 
^ , ing) forming on the "background" of these pseudogap fluctuations. Microscopic 

derivation of Ginzburg - Landau expansion allows calculations of critical temper- 
ature Tc and other basic characteristics of a superconductor, depending on the 
parameters of the pseudogap. We also analyze the role of "normal" (nonmagnetic) 
^ I impurity scattering. It is shown that our simplified model allows semiquantitative 

■ modelling of the typical phase diagram of superconducting cuprates^ . 



^Extended version of the talk given by the author on the seminar of I.E.Tamm Theoretical Department 
of P.N.Lebedev Physical Institute, Russian Academy of Sciences, Moscow, October 7, 2003. 
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Figure 1: Typical phase diagram of copper oxide high - temperature superconductor. 
Different phases and anomalous regions are shown schematically. Carrier concentration, 
corresponding to a maximum of transition temperature is called "optimal". Systems 
with lower concentrations are called "underdoped" , those with larger — "overdoped" . 

1 Basic problems of the physics of high — tempera- 
ture superconductors. 

We shall start with brief review of the present day situation in the physics of high - 
temperature superconductors, which may be useful for the reader not involved in this 
field. 

1.1 What is really KNOWN about copper oxides: 

High - temperature superconducting (HTSC) copper oxides are intensively studied for 
more than 15 years now. In these years, considerable progress has been achieved in our 
understanding of the nature and basic physical properties of these systems, despite their 
complicated phase diagram, containing almost all the main phenomena studied by solid 
state physics (Fig. H}. 

It is apparent that the number of facts are now well established and are not due to 
any revision in the future ^ If we talk about superconducting state, the list of "solved" 
problems may be formulated as follows: 

• Nature of superconductivity ^ Cooper pairing. 

^Surely, any such list is subjective enough and is obviously based on some prejudices and "tastes" of 
the author. It is well known, that in HTSC physics quite opposite views are often "peacefully coexistent" 
with each other. However, such list may hopefully be of some interest. 
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It is definitely established that in these system we are dealing with 

1. d - wave pairing, 

and in the spectrum of elementary excitations we observe 

2. Energy gap A cos 20, 

where the polar angle determines the direction of electronic momentum in 
two - dimensional inverse space, corresponding to Cu02 plane. Thus, energy 
gap becomes zero and changes sign on the diagonals of two - dimensional 
Brillouin zone. 

3. The size of pairs is relatively small: ~ 5 — 10a, 

where a is the lattice constant. This means that high - temperature supercon- 
ductors belong to a crossover region between "large" pairs of BCS theory and 
"compact" Bosons picture of very strongly coupled electrons. 

Superconducting state always appears as 

4. Second order phase transition 

with more or less usual thermodynamics. 

All these facts were established during approximately first five years of studies of 
HTSC copper oxides. We do not give any references to original works here as the 
list will be to long, and quote only two review articles 

Situation becomes much more complicated when we consider the properties of the 
normal state. Almost all anomalous properties of these unusual systems appear first 
of all in the normal state. However, here we also can list some definitely established 
facts: 

• Existence of the Fermi surface. 

In this sense these systems are clearly metals, though rather unusual ("bad"). Most 
of the data on Fermi surfaces were obtained from ARPES and the progress 

here in recent years is quite spectacular due to the great increase (more than order 
of magnitude!) in resolution both in energy and momentum. As an illustration 
of these advances in Fig. |21 we show experimentally determined Fermi surface of 
La2-xSrxCuOi with x = 0.063 0. Note the existence of characteristic "fiat" parts 
of this Fermi surface. This form of the Fermi surface is rather typical for the majority 
of HTSC cuprates in superconducting region of the phase diagram. 

• Metal — Insulator transition, 

which takes place with the change of chemical composition (as the number of carriers 
drops). Stochiometric La2Cu04 is an antiferromagnetic insulator (apparently of 
Mott type) with Neel temperature of the order of 400K. There exists well defined 
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Figure 2: Experimentally determined Fermi surface of La2~xSrxCu04 with x = 0.063. 
At the insert - schematic intensity of neutron scattering with four incommensurate peaks 
close to the scattering vector (tt, tt). Dashed lines show the borders of antiferromagnetic 
Brillouin zone, which appear after the establishment of antiferromagnetic long - range 
order (period doubling). Points, where these borders intersect the Fermi surface are 
called "hot" spots. 
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optical (energy) gap, and antiferromagnetism is due to the ordering of localized spins 
on Cu ions and is well described by two - dimensional Heisenberg model. This is 
more or less typical also for other HTSC oxides, which thus belong to a wide class 
of strongly correlated electronic systems. At the same time, this dielectric state is 
rapidly destroyed by introduction of few percents of doping impurities. Note that 
the Fermi surface shown in Fig. |2l is observed, in fact, in the system which is on 
the edge of metal - insulator transition, which takes place at x ~ 0.05. 

• Strong anisotropy of all properties (quasi - two - dimensional nature!). 

Current carriers (in most cases — holes) are propagating more or less freely along 
Cu02 planes (a and h directions in orthorombic crystal), while transverse motion 
in orthogonal direction (c axis) is strongly suppressed. Conductivity anisotropy 
is usually of the order of 10^ — 10^. This fact is seriously damaging from the 
point of view of practical applications. At the same time it is still unclear whether 
two - dimensionality is a necessary condition for realization of high - temperature 
supeconductivity. 

1.2 What is still UNKNOWN about copper oxides: 

Let us now list the main unsolved problems, remaining in the center of rather sharp 
discussions with participants sometimes not hearing each other. On the first place stays 
of course the 

• Mechanism of Cooper pairing. 

Most researchers (not all!) do not have any doubt, that in HTSC - oxides we are 
dealing with Cooper pairing within more or less standard BCS "scenario". However, 
the question is — what kind of interaction leads to pair formation? A number of 
microscopic mechanisms are under consideration: 

1 . Electron - phonon j^l E] . The main unsolved problem here is an explanation 
of d - wave symmetry of Cooper pairing. Usually, electron - phonon coupling 
leads to s - wave pairing. Recently there was some progress in possible ways 
to solve this problem [H]. 

2. Spin - fluctuation CHI CIl IIS • 

Historically, within this approach the d - wave nature of pairing in oxides 
was predicted. Further experimental confirmation, as well as the possibility 
of semiquantitative description of many properties of these systems made this 
approach probably most favorable (from my point of view!) among many other 
mechanisms. 
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3. Exchange - RVB, S0{5), ...? 

These models ITl] , as well as a number of more "exotic" , were formulated 
in early days of HTSC research. Enormous intellectual resources of leading 
theorists were and are spent on their development. However (again it is only 
my personal view), the completeness and "stability" of results here is incom- 
parably less than in more traditional approaches, while real connection with 
experiments is relatively slight, so that all these models still remain a kind 
of "brain gymnastics". Obviously, I do not deny usefulness of these models 
from purely theoretical point of view. For example, the discussion of symme- 
try connections of antiferromagnetism and superconductivity ^3] is of great 
interest. At the same time, the most natural and useful among phenomenolog- 
ical approaches remains the description of superconductivity in oxides based 
on anisotropic version of the standard Ginzburg - Landau theory ^. 

At this time it is unclear, which of these mechanisms is realized or is dominant in 
real HTSC - oxides. And the reason for this lies in the fact, that the properties of 
superconducting state are relatively independent of microscopic mechanism of pair- 
ing. It is rather difficult to propose some "crucial" experiment, which will definitely 
confirm one of these mechanisms. In some sense, situation here is somehow similar 
to those existing, actually for decades, in the theory of magnetism. Since the mid- 
dle of the thirties, it is well known, that the nature of magnetism is connected with 
exchange interaction. There exist plenty microscopic mechanisms of exchange inter- 
action between spins (e.g. direct exchange, super exchange, s — d interaction, RKKY 
etc.). However, it is not always possible to tell, which of these mechanisms acts in 
some real magnetic system. Classical example here is the problem of magnetism of 
iron! 

• Nature of the normal state. 

Here any consensus among researchers is almost absent. Usually, recognizing the 
metallic nature of these systems, theorist discuss several alternative possibilities: 

1. Fermi - liquid (Landau). 

2. "Marginal" or "bad" Fermi - liquid [H]. 

3. Luttinger liquid |13j. 

All the talking on the possible breaking of Fermi - liquid behavior (absence of 
"well defined" quasiparticles) originates, from theoretical point of view, from 
strong electronic correlations and quasi - two - dimensionality of electronic 
properties of these systems, while from experimental side this is mainly due 

■^This fact is of special significance as this seminar has taken place approximately an hour after the 
news came on the award of 2003 Nobel Prize in physics to V.L. Ginzburg. 
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to practical impossibility of performing studies of electronic properties of the 
normal state at low enough temperatures (when these properties are in fact 
"shunted" by superconductivity, which is impossible to suppress without sig- 
nificant changes of the system under study). One must remember, that Fermi - 
liquid behavior, by definition, appears only as limiting property at low enough 
temperature. Thus, it is not surprising at all that it is "absent" in experiments 
performed at temperatures of the order of lO^K! At the same time, as we 
shall see below, recently there were an important experimental developments, 
possibly clarifying the whole problem. 

4. Disorder and local inhomogeneities. 

Practically all high - temperature superconductors are internally disordered 
due to chemical composition (presence of doping impurity). Thus, any under- 
standing of their nature is impossible without the detailed studies of the role 
of an internal disorder in the formation of electronic properties of strongly cor- 
related systems with low dimensionality. In particular, localization effects are 
quite important in these systems, and studies of these has already a long his- 
tory ^TC]. Situation with the role of internal disorder has complicated recently 
with the arrival of the new data, obtained in by scanning tunneling microscopy 
(STM), which clearly demonstrated inhomogeneous nature of the local density 
of states and superconducting energy gap on microscopic scale, even in prac- 
tically ideal single - crystals of copper oxides. Of many recent papers devoted 
to these studies, we quote only two [T7[ HH], where further references can be 
found. These results significantly modify previous ideas on microscopic phase 
separation, "stripes" etc. On the other hand, the presence of such inhomo- 
geneities makes these systems a kind of "nightmare" for theorists, though the 
picture of inhomogeneous superconductivity due to fluctuations in the local 
density of states was analyzed rather long time ago [TBI IT^. 

5. Pseudogap. 

Now it is clear that most of the anomalies of the normal state of copper oxides is 
related to the formation of the so called pseudogap state. This state is realized 
in a wide region of the phase diagram as shown in Fig. ^ corresponding mainly 
to "underdoped" compositions and temperatures T < T*. It is important to 
stress, that the line T* on Fig. ^ is rather approximate and apparently do 
not correspond to any phase transition and signify only a kind of crossover 
to the region of well developed pseudogap anomahes^. In short, numerous 
experiments [201 1^ show that in the region of T < T* the systems, which are 
in the normal (non superconducting) state, do possess some kind of the gap - 

^Sometimes the notion of "high energy pseudogap" is introduced and defined by T*, while the "low 
energy pseudogap" is defined by another crossover line in the region of T < T* closer to superconducting 
"dome" , as shown in Fig. 
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like feature in the energy spectrum. However, this is not a real gap, but some 
kind of precursor of its appearance in the spectrum (which is the reason for 
the term "pseudogap" , which appeared first in qualitative theory of amorphous 
and liquid semiconductors |22j.). Some examples of the relevant experiments 
will be given below. 

Crudely speaking, the gap in the spectrum can be either of superconducting or 
insulating nature. Accordingly, there exist two possible theoretical "scenarios" 
to explain pseudogap anomalies in copper oxides. The first one anticipates 
formation of Cooper pairs already at temperatures higher, than the temper- 
ature of superconducting transition, with phase coherence appearing only in 
superconducting state at T < Tc. The second assumes, that the origin of the 
pseudogap state is due to fluctuations of some kind of short range order of 
"dielectric" type, developing in the underdoped region. Most popular here is 
the picture of antiferromagnetic (AFM) fluctuations, though fluctuating charge 
density waves (CDW) cannot be excluded, as well as structural deformations 
or phase separation at microscopic scales. In my opinion, most of the recent 
experiments provide evidence for this second scenario. So, in further discussion 
we shall deal only with this type of models, mainly speaking about antiferro- 
magnetic fluctuations. 

It is common view, that the final understanding of the nature of high - temperature 
superconductivity is impossible without clarification of the nature of the normal state. 
Generally speaking, we have to explain all characteristic features of the phase diagram, 
shown in Fig. ^ This is the main task of the theory and it is still rather far from the 
complete solution. In the following we shall concentrate on the discussion of some simple 
models of the pseudogap state and attempts to illustrate possible ways to find a solution 
of this main problem. 

2 Basic experimental facts on the pseudogap behav- 
ior in high — temperature superconductors. 

Phase diagram of Fig. shows that, depending on concentration of current carriers in 
highly conducting Cu02 plane, a number of phases and regions with anomalous proper- 
ties can be observed in copper oxides. For small concentrations, all the known HTSC - 
cuprates are antiferromagnetic insulators. With the growth of carrier concentration Neel 
temperature T/v rapidly drops from the values of the order of hundreds of vanishing 
at concentration of holes p less than or of the order of 0.05, and system becomes metallic. 
As concentration of holes is increasing further, the system becomes superconducting, and 
temperature of superconducting transition grows with the growth of carrier concentra- 
tion, passing through characteristic maximum at po ~ 0.15 — 0.17 (optimal doping), and 
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then drops and vanish at p ~ 0.25 — 0.30, though metalhc behavior in this (overdoped) 
region remains. In the overdoped region with p > Po metalhc properties are more or less 
traditional (Fermi - liquid), while for p < po the system is a kind of anomalous metal, not 
described (in the opinion of majority of authors) by the usual theory of Fermi - liquid. 

Anomalies of physical properties attributed to the pseudogap formation are observed in 
metallic phase for p < p^ and temperatures T < T*, where T* drops from the values of the 
order of T/v at p ~ 0.05 vanishing at some "critical" concentration of carriers Pc, slightly 
greater thanpo, so that T* - line in Figgis actually continuing under the superconducting 
"dome". According to the data presented in Ref. |23j T* vanishes at p = pc ^ 0.19. In the 
opinion of some other authors (proponents of superconducting scenario for the pseudogap) 
T* - line just passes to - line of superconducting transition somewhere close to the 
optimal doping Pq. In my opinion, most recent data, apparently, confirm the first variant 
of the phase diagram (more details can be found in Ref. pBj ) . 

Pseudogap anomalies, in general, are interpreted as due to suppression (in this region) 
of the density of single - particle excitations close to the Fermi level, which corresponds 
to the general concept of the pseudogap The value of T* then determines charac- 
teristic scale of the observed anomalies and is proportional to effective energy width of 
the pseudogap. Let us now consider typical experimental manifestations of pseudogap 
behavior. 

Consider first experimental data on electronic specific heat of cuprates. In metals this 
contribution is usually written as: C = 'y(T)T, so that in the normal state (T > Tc) 
7 ~ N{0), where A^(0) - is the density of states at the Fermi level. At T = Tc we have 
the well known anomaly due to second order phase transition, so that 7(T) demonstrates 
characteristic peak (discontinuity). As a typical example, in Fig. El we show experimental 
data for Yo,sCao,2Ba2Cu30j-s for different values of 6 We can see, that in optimally 
doped and overdoped samples 7(T) is practically constant for all T > Tc, while for 
underdoped samples considerable drop of 7(T) appears for T < 150 — 200K. This is a 
direct evidence of the appropriate drop of the density of states at the Fermi level due to 
pseudogap formation for T < T*. 

Note also that the value of specific heat discontinuity at superconducting Tc is signif- 
icantly suppressed as we move to the underdoped region. More detailed analysis shows 
j2S| that the drop of this discontinuity A'jc actually starts at some "critical" carrier con- 
centration Pc ~ 0.19, which is connected with the "opening" of the pseudogap. 

Pseudogap formation in the density of states is clearly seen also in experiments on 
single particle tunneling. Thus, in highly cited Ref. |2S] tunneling experiments were per- 
formed on single crystals of Bi2Sr2CaCu20s+s {Bi — 2212) with different oxygen content. 
For underdoped samples pseudogap formation in the density of states was clearly observed 
at temperatures significantly higher than superconducting Tc. This pseudogap smoothly 
transformed into superconducting gap at T < Tc, which is often seen as an evidence of its 
superconducting nature. However, in Refs. [231211; where tunneling experiments were 
performed on the same system, it was directly shown, that superconducting gap exists 
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Figure 3: Coefficient 7 of electronic specific heat in overdoped (a) and underdoped (b) 
Yo,8Cao,2Ba2Cus07_s ■ 

on the "background" of the wider pseudogap and vanishes at T = Tc, while pseudogap 
persists at higher temperatures. 

Pseudogap also shows itself in transport properties of HTSC - systems in normal state, 
Knight shift and NMR relaxation. In particular, changes from the usual for optimally 
doped linear temperature dependence of resistivity for underdoped samples and T < T* 
are often attributed to its formation. Also the value of Knight shift in such samples for 
T < T* becomes temperature dependent and drops as temperature lowers. Similar behav- 
ior is observed in underdoped samples for (TTi)~^, where Ti - is NMR relaxation time. 
Let us remind, that in usual metals Knight shift is just proportional to the density of states 
at the Fermi level A^(0), while (TTi)~^ ~ A^^(O) (Korringa behavior), and resistivity p is 
proportional to the scattering rate (inverse mean free time) 7 ~ N{0). Thus, significant 
lowering of these characteristics is naturally attributed to the drop in the density of states 
A^(0) at the Fermi level. Note that these arguments are, of course, oversimplified, partic- 
ularly when we are dealing with temperature dependences. E.g. in case of resistivity this 
dependence is determined by inelastic scattering and physics of these processes in HTSC 
is still unclear. Thus, the decrease in the density of states (due to partial dielectrization 
of the spectrum) can also lead, in fact, to the growth of resistivity. 

Pseudogap in underdoped cuprates is also observed in experiments on optical con- 
ductivity, both for electric field polarization along highly - conducting CUO2 plane and 
also along orthogonal direction (c - axis). These experiments are reviewed in detail in 
Ref. [23 . As a typical example in Fig. |3]we show the data of Ref. j2H] on optical con- 
ductivity in the CUO2 plane for different compositions of electronically conducting oxide 
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Nd2-xCexCuOi- Characteristic feature here is the appearance (most clearly seen for un- 
derdoped samples with x = 0.05 and x = 0.10), of the shallow minimum at frequencies of 
the order of 0.25 eV and of smooth maximum due to absorption through unusually wide 
pseudogap around u ~ 0.5 eV. The nature of an additional absorption maximum seen at 
a; ~ 0.1 eV can, apparently, bw attributed to localization (cf. below). 

In underdoped cuprates with hole conductivity, optical conductivity is usually char- 
acterized by a narrow "Drude like" absorption peak at small frequencies, followed by a 
shallow minimum and smooth maximum due to pseudogap absorption through pseudo- 
gap of the order of O.leV j20| I21j. Additional peak due to localization is observed rarely, 
usually after the introduction of an additional disorder [211 EHl HI] • 

Most spectacular effects due to pseudogap formation are seen in experiments on angle 
resolved photoemission (ARPES). ARPES intensity (energy and momentum distribution 
of photoelectrons) is determined by 

likcu) = /o(k)/(^)A(ka;) (1) 

where k - is the momentum in the Brillouin zone, u) - energy of initial state, measured 
with respect to the Fermi level (chemical potential)^, /o(k) includes some kinematic factors 
and the square of the matrix elements of electron - photon interaction and in crude 
approximation is considered to be some constant, 

A(k,u) = --ImG(k,Lj + i6) (2) 

TT 

where G{k, uj) - is Green's function of an electron, determines the spectral density. The 
presence of Fermi distribution /(tu) = [exp{uj/T) + 1]^^ reflects the fact, that only occu- 
pied states can produce photoemission. Thus, in such a crude approximation it can be 
said, that ARPES experiments just measure the product f{u)A{kuj), and we get direct 
information on the spectral properties of single - particle excitations. 

Consider qualitative changes in single - electron spectral density Q due to pseudogap 
formation. In the standard Fermi - liquid theory, single - electron Green' function can 
be written as: 

G{UJ, k) = : h Gincoh (3) 

- 4k - «7k 

where = ~ ^ is quasiparticle energy with respect to the Fermi level (chemical 
potential) /i, 7k - quasiparticle damping. The residue in the pole < Zk < 1, Gincoh ^ 
is some non singular contribution due to many particle excitations. Then the spectral 
density is: 

Aiu, k) = -Z^- \ , + ... (4) 

TT (a;-^k)^ + 7k 

^In real experiments lo is measured with respect to the Fermi level of some good metal, e.g. Pt or Ag, 
placed in electric contact with a sample. 
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Figure 4: Real part of optical conductivity in the ab - plane of Nd2-xCexCu04 for 
different temperatures and compositions x. 
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where dots denote more or less smooth contribution from Gincoh, while quasiparticle spec- 
trum determines a narrow (if damping 7^ is small in comparison to ^k) Lorentzian peak. 
In the usual Fermi - liquid 7 ~ u)"^ ~ |k — kiT'p and quasiparticles are well defined in 
some (narrow enough) vicinity of the Fermi surface. In the model of "marginal" Fermi - 
liquid 7 ~ cij ~ |k — hpl and quasiparticles are "marginally" defined due to 7 ~ ^k- In 
the presence of static scattering (e.g. due to impurities) a constant (independent of uj) 
contribution 70 appears in damping. 

Qualitative form of the spectral density in the Fermi - liquid picture is illustrated in 
Fig. Eta). 

If some long - range order (e.g. of SDW(AFM) or CDW type) appears in the system, 
an energy gap Wk (of dielectric nature) opens in the spectrum of elementary excitations 
(k - dependence stresses the possibility of gap opening only on the part of the Fermi 
surface), and single - particle Green's function acquires the form of Gorkov's function 
(where we also add some damping Fk): 

^(^,1^) = ^ + (5) 

u; - Ek + iFk + Ek - "^rk 

where the excitation spectrum is now: 

= + (6) 
and we introduced Bogoliubov's coefficients: 



< = -J^ + ^ (7) 




Then the spectral density is: 



4/ 1 ^ Tk fk - K /p,N 



where now appear two peaks, narrow if Fk is small enough, corresponding to "Bogoli- 
ubov's" quasiparticles. 

If there is no long - range order, but only strong scattering by fluctuations of short - 
range order, characterized by some correlation length ^, is present in our system, it is easy 
to imagine, that spectral density possesses (in the same region of momentum space and 
energy) some precursor structure, in the form of characteristic "double - hump" structure, 
as it is shown qualitatively in Fig. El The widths of these maxima are naturally determined 
by parameter Vp/C,, i.e. inverse time of flight of an electron through the region of the size 
of ^, where "dielectric" ordering is effectively conserved. Below we shall see that rigorous 
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Figure 5: Qualitative evolution of spectral density, (a) - normal metal (Fermi - liquid), 
= - at the Fermi surface, (b) - two narrow peaks, corresponding to "Bogoliubov's" 

quasiparticles in a system with dielectric gap (in case of long - range order of CDW or 

SDW type). Smooth maxima - system without long - range order (pseudogap behavior), 
= - at the Fermi surface, (c) - same as (b), but for ,^k > 0, i.e. above the Fermi 

surface. Note characteristic asymmetry of maxima in this case. 
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Figure 6: ARPES spectra in three different points a, b, c on the Fermi surface of Bi — 2212 
for underdoped sample with Tc = 85K. Thin curves — spectrum of the reference sample 
of Pt. 

analysis leads just to these results. In this sense, the further theoretical discussion will 
be devoted to justification of this qualitative picture. 

In Fig. iniwe show ARPES data for Bi2Sr2CaCu20s |22j, obtained in three different 
points on the Fermi surface for different temperatures. The presence of the gap (pseudo- 
gap) manifests itself by the shift (to the left) of the leading edge of energy distribution of 
photoelectrons in comparison to the reference spectrum of a good metal (Pt). It is seen 
that this gap is closed at different temperatures for different values of k, and the gap width 
diminishes as we move from (0, 0) — (0, vr) direction in the Brillouin zone. Pseudogap is 
completely absent in the direction of zone diagonal (0,0) — (vr, vr). At low temperatures 
this is in complete accordance with the picture of - wave pairing, which is confirmed 
in cuprates by many experiments QUI]- Important thing, however, is that this "gap" in 
ARPES data is observed also at temperatures significantly higher than the temperature 
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of superconducting transition T^.. 

In Fig. [7|we show angular dependence of the gap in the Brillouin zone and temperature 
dependence of its maximal value, obtained from ARPES for several samples of Bi — 
2212 with different compositions. It is seen that the general d - wave like symmetry 
is conserved and the gap in optimally doped sample vanishes practically at T = Tc, 
while for underdoped samples we observe typical "tails" in the temperature dependence 
of the gap for T > T^. Qualitatively we may say that the formation of an anisotropic 
pseudogap at T > Tc, which smoothly transforms into superconducting gap for T < Tc, 
leads to "destruction" parts of the Fermi surface of underdoped samples, close to the 
point (0, vr) (and symmetrical to it), already for T < T*, and the sizes of these parts grow 
as temperature lowers [22]. 

Of central interest is, of course, the evolution of the spectral density A{\ip,uj) at the 
Fermi surface. Under rather weak assumptions it can be directly determined from ARPES 
data [32] ■ Assuming electron - hole symmetry (always valid close enough to the Fermi 
surface, in reality for |a;| less than some tenths of meV) we have A(^Fi^) = Aikp, —u), 
so that taking into account f{—uj) = 1 — f{—uj), from and for k = k^T^ we obtain 
I{uj) +/(— cu) = yl(k^,a;). Thus, the spectral density at the Fermi surface can be directly 
determined using symmetrized experimental data I{oj) + As an example of such 

analysis, in Fig. |H] we show the dtat of Ref. [31] for underdoped sample of Bi — 2212 
with Tc = 83K and overdoped with Tc = 82/^ at different temperatures. It is seen that 
pseudogap existence clearly manifests itself in characteristic "double - humps" structure of 
spectral density, which appears (in an underdoped system) for temperatures significantly 
higher than T^.. We see that these data completely correspond to the expected qualitative 
form of spectral density in the pseudogap state. 

Let us stress once again, that well defined quasiparticles correspond to narrow enough 
peak in the spectral density A{\iFi^) at = 0. Such behavior, until recently, was in 
practically never observed in copper oxides. However, it was discovered some time ago 
that in superconducting phase, at T <C Tc, there exists sharp enough peak of the spectral 
density, corresponding to well defined quasiparticles, in the vicinity of an intersection of 
the Fermi surface with diagonal of the Brillouin zone (direction (0,0) — (7r,7r)), where 
superconducting gap vanishes ^]. At the same time, close to the point (0, tt) the Fermi 
surface remains "destroyed" by both superconducting gap and the pseudogap. The studies 
of these "nodal" quasiparticles id of great importance and lead to some clarification of the 
problem of Fermi - liquid behavior. Most recent data show, that quasiparticle peak in 
diagonal direction persists also at temperatures much higher than Tc. This is clearly seen 
from the data of Ref. [SH] , shown in Fig. [21 where we can see the evolution of this peak as 
we move along the Fermi surface. It is seen that quasiparticle behavior is valid everywhere 
for overdoped samples and only in the vicinity of diagonal for underdoped (and optimal). 
Quasiparticle peak in diagonal direction persists even for strongly underdoped samples, 
which are on the edge of metal - insulator transition [3] . From data shown in Fig. ^3B] , 
we can also see significant anisotropy of static (or more precisely quasistatic within limits 



17 




O 50 lOO 150 200 250 300 



T(K) 

Figure 7: Angular and temperature dependences of energy gap in Bi — 2212, obtained 
from ARPES data for samples with = 87 K (nearly optimally doped), = 83K and 
Tc — lOK (underdoped) : (a) - the value of ARPES gap, measured at different points of the 
Fermi surface (shown at the insert), with positions determined by polar angle, measured 
from the direction FM, d - wave symmetry of the gap is obvious, (b) - temperature 
dependence of the maximal gap, measured close to point M. 
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Figure 8: Symmetrized ARPES spectra for overdoped sample of Bi — 2212 with Tc = 827^^ 
(a) and underdoped sample with Tc = 83K (b) at the point of intersection of the Fermi 
surface with the border of Brillouin zone (0, n) — (vr, vr). 

of ARPES resolution) scattering, growing as we move to the vicinity of (0, tt), which can 
obviously be related to pseudogap formation. 

Finally, in Fig. ^2 we show spectacular results of Ref . jSTj , obtained from ARPES 
measurements with significantly improved resolution, from which, even by the "naked 
eye", we can see that effective damping of "nodal" quasiparticles is in general (except 
possibly underdoped samples) more or less quadratic in energy, as we can expect for a 
standard Fermi - liquid. Only close to optimal doping and well into underdoped region 
we can see some additional contribution to scattering, most probably of magnetic nature. 
These newest results provide quite new information on much discussed problem of Fermi 
- liquid behavior in HTSC - oxides. Apparently we observe typical Fermi - liquid be- 
havior in overdoped samples, but it is "destroyed" as we move to optimal doping and 
into underdoped region, though only in parts of the momentum space (in the vicinity of 
the points like (0,7r)), where pseudogap appears, leading to additional strong scattering 
(damping). 

The growth of pseudogap anomalies as we move to the region around (0, vr), where also 
becomes maximal the amplitude of superconducting d - wave gap, is often interpreted as 
an evidence of d - wave symmetry of the pseudogap and of its superconducting nature, 
smoothly transforming into the real gap for T < Tc. However, there is lot of evidence that 
pseudogap actually competes with superconductivity and is, most probably, of dielectric 
nature. Detailed discussion of this evidence can be found e.g. in Refs. j^Il 123 • Here 
we only limit ourselves to the most general argument — pseudogap anomalies in HTSC 
grow as we move deep into underdoped region, where superconductivity just vanish. It is 
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Figure 9: ARPES spectra, measured at different point at the Fermi surface, moving from 
the diagonal direction in Brillouin zone to (vr, 0): (a) - points, where data were taken; (b) 
- comparison of data in "antinodal" point A (vicinity of (tt, 0)) in optimally doped and 
overdoped Bi22Ql; (c) - data for optimally doped Bi2212 with Tc — 90K, obtained at 
T = UOK; (d) - similar data for overdoped Bi2201 with = 0, obtained at T = UOK. 
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Figure 10: The fit of data on ARPES spectrum width to "marginal" dependence of the 
imaginary part of self - energy: 7mE — a + bcu. Clearly seen is anisotropy of static 
(quasistatic) scattering in optimally doped Bi2212, while dynamic (inelastic) scattering 
remains isotropic (c). Behavior of the velocity on the Fermi surface as shown on (d) and 
(e) (overdoped Bi2201). 
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Figure 11: Experimental data on the width of ARPES peaks in Bi{Pb)2212 and Bi2212 for 
different temperatures, as doping level changes from underdoped sample with Tc = 76K 
(a) to overdoped with = 73K (d). 



difficult to imagine, how this type of behavior can be due to precursor Cooper pairing at 
T>Te. 

Of course, there are experimental data, giving direct evidence for "dielectric" nature 
of the pseudogap. In Fig. we present ARPES data on the Fermi surface of elec- 
tronic superconductor Ndi s^Ceo i^CuO^ jSH], which clearly show, that the "destruction" 
of the Fermi surface due to pseudogap formation takes place in the region around "hot 
spots", appearing at intersections of the Fermi surface with borders of the "future" an- 
tiferromagnetic Brillouin zone, which would have appeared after the establishment of 
antiferromagnetic long - range order. Real gap in this case is obviously of the usual 
"band" or "insulating" nature, nothing to do with Cooper pairing of d - wave symmetry. 
Most ARPES experiments in HTSC - oxides are performed on systems with hole - like 
conductivity, where the distance (in momentum space) between sides of the Fermi surface 
in the vicinity of (0, tt) is just smaller, than in NdCeCuO . Also smaller in these systems is 
characteristic energy scale (width) of the pseudogap, as was mentioned previously during 
the discussion of optical data shown in Fig. 0] Accordingly, the resolution of ARPES is 
apparently insufficient to resolve separate "hot spots" , which are too close to each other in 
the vicinity of (0, vr). In my opinion, the results of Ref. jSHI practically solve the problem 
in favor of "dielectric" scenario of pseudogap formation in cuprates. 

In conclusion, we must stress that in different experiments discussed above, charac- 
teristic temperature T*, defining crossover into the pseudogap state can somehow change, 
depending on the property which is being studied. However, in all cases there is some 
systematic dependence of T* on the doping level and this temperature vanishes at some 
concentration of carriers slightly higher than optimal. In Fig. ^Jwe show a compendium 
of data (derived from a number of different experiments) on the energy width of the 
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Figure 12: Fermi surface of electronic superconductor Ndi s5Ceo,i5Cu04 in one quadrant 
of the Brillouin zone, obtained by direct integration of ARPES spectra. Dashed hne 
denotes the border of antiferromagnetic Brillouin zone. 




Figure 13: Dependence of the energy width of the pseudogap Eg in YBCO on hole 
concentration determined from different experiments. 
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pseudogap Eg for YBCO as a function of hole concentration [53] (in this reference it was 
assumed, rather arbitrarily, that Eg = 2.5T*). It is seen that pseudogap vanishes at the 
critical concentration pc ~ 0.19, slightly higher than optimal concentration of carriers, 
which, by the way, is an additional argument in favor of non superconducting nature of 
the pseudogap in cuprates. 

More detailed review of miscellaneous experiments on the pseudogap formation in high 
temperature superconductors can be found in Refs. [201 122] ■ 

3 Theoretical considerations — variants of simplified 
model. 

As we mentioned above, there are two alternative scenarios to explain pseudogap anoma- 
lies in HTSC - systems. The first one is based on the model of Cooper pair formation 
already above the temperature of superconducting transition (precursor pairing), while the 
second assumes, that the origin of the pseudogap state is due to scattering by fluctuations 
of short - range order of "insulating" type (e.g. antiferromagnetic (SDW) or charge den- 
sity wave (CDW)), existing in underdoped region of copper oxides. This second scenario 
seems to be more attractive, both due to a number of experimental evidences and due to a 
simple fact, that all pseudogap anomalies become stronger, as carrier concentration drops 
and the system moves farther away from optimal concentration for superconductivity 
towards dielectric (antiferromagnetic) phase. 

Consider typical Fermi surface of electrons moving in the Cu02 plane, shown in Fig. 
121 If we neglect fine details, the observed (e.g. in ARPES) Fermi surface topology (and 
also the spectrum of elementary excitations) in Cu02 plane, in the first approximation 
are well enough described by the usual tight - binding model: 

= — 2t(cos k^a + cos kya) — At' cos k^a cos kya (10) 

where t ^ 0.25eV is the nearest neighbor transfer integral, while t' is the transfer in- 
tegral between second - nearest neighbors, which can change between t' ^ — 0.45t for 
YBa2Cu30T_s and f ^ — 0.25t for La2-xSrxCuOi, a is the square lattice constant. 

Phase transition to antiferromagnetic state induces lattice period doubling and leads 
to the appearance of "antiferromagnetic" Brillouin zone in inverse space as shown in Fig. 
[Til If the spectrum of carriers is given by ()10p with t' = and we consider the half - filled 
case, Fermi surface becomes just a square conciding with the borders of antiferromagnetic 
zone and we have a complete "nesting" — flat parts of the Fermi surface match each 
other after the translation by vector of antiferromagnetic ordering Q = (±7r/a, zbvr/a). 
In this case and for T = the electronic spectrum is unstable, energy gap appears 
everywhere on the Fermi surface and the system becomes insulator, due to the formation 
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of antiferromagnetic spin density wave (SDW)^. This picture corresponds to one of the 
popular schemes to explain antiferromagnetism in cuprates, see e.g. Ref. [SS] and a 
review paper jlHI) though, as was noted above, experimental situation is well described 
by simple Heisenberg model with localized spins. In the case of the Fermi surface shown in 
Fig^^the appearance of antiferromagnetic long - range order, in accordance with general 
rules of the band theory, leads to the appearance of discontinuities of isoenergetic surfaces 
(e.g. Fermi surface) at crossing points with borders of new (magnetic) Brillouin zone due 
to gap opening at points connected by vector Q. Note similarity of this picture with 
experimental data shown in Fig. ^1 

In the region of cuprate phase diagram of interest to us antiferromagnetic long - range 
order is absent, however, a number of experiments support the existence (everywhere 
below the T* - line) of well developed fluctuations of antiferromagnetic short - range 
order which scatter electrons with characteristic momentum transfer of the order of Q. 
In principle, it is not very important to consider AFM(SDW) fluctuations, similar effects 
may be due CDW fluctuations. There is no doubt, that such scattering processes exist 
(and dominate!) in cuprates in the whole pseudogap region. To convince yourself, just 
look at the data shown in Fig. EI 

For concreteness consider, however, the model of "nearly antiferromagnetic" Fermi 
- liquid 02]; where the effective interaction of electrons with spin fluctuations is 
described by dynamic spin susceptibihty Xq(a;), the form of which was determined by 
fitting to NMR experiments Ej : 



where g is coupling constant, ^ - correlation length of spin fluctuations, Q = (±7r/a, ±n/a) 
is vector of antiferromagnetic ordering in dielectric phase, Ugf ^ characteristic frequency 
of spin fluctuations. 

Dynamical spin susceptibility x^{^) is peaked around wave vectors (±7r/a, ±7r/a), 
and this leads to appearance of "two types" of quasiparticles — "hot" one, with momenta 
in the vicinity of "hot spots" on the Fermi surface fFiglT^ and energies satisfying the 
inequality {vp - velocity at the Fermi surface): 



and "cold" one with momenta close to the parts of the Fermi surface surrounding diagonals 
of Brillouin zone \px\ = \py\ and not satisfying ()12|). This terminology is connected with 
strong scattering of quasiparticles in the vicinity of "hot spots" with momentum transfer 
of the order of Q due to interaction with spin fluctuations (jlip . while for quasiparticles 
with momenta far from "hot spots" this interaction is weak enough. In the following 

^Analogous dielectrization is realized also in the case of the formation of the similar charge density 
wave (CDW). 



Veff{q,Uj) =/Xq(^) 




(11) 




(12) 
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Figure 14: Fermi surface in the Brillouin zone and "hot spots" modeL Magnetic zone 
appears in the presence of antiferromagnetic long - range order. "Hot spots" correspond 
to intersections of its borders with Fermi surface and are connected by the scattering 
vector of the order of Q = (-, - ). 
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we shall call this "hot spots" model. Correlation length of fluctuations of short - range 
antiferromagnetic order ^, described by (llip . is an important parameter of this theory. 
Note that in real HTSC - systems this length is not very large, usually 2a < C, < 8a 



Characteristic frequency of spin fluctuations cUgf, depending on the compound and 
doping level, is usually in the limits of 10 — lOOK HEj; so that in most part of the 
pseudogap region on the phase diagram we have 27rT ^ oUsf and actually can neglect spin 
dynamics, limiting ourselves to quasistatic approximation: 

where W is an effective parameter with dimensions of energy, which in the model of AFM 
fluctuations can be written as [TF] : 

W = = g'< {n.; - n,,r > (14) 

where g is interaction constant of electrons and spin fluctuations, < > is the aver- 
age square of spin on a lattice site, nj|, - operators of a number of electrons on a 
given site with appropriate spin direction. In this approximation, dynamic fleld of spin 
fluctuations is just replaced by the static Gaussian random field of "quenched"^ spins, 
antiferromagnetically correlated on lengths of the order of ^. 

It is clear that in the framework of our semiphenomenological approach, both corre- 
lation length ^ and parameter W are to be considered as some functions of carrier con- 
centration (and temperature) to be determined from the experiment. In particular, W 
determines the effective width of the pseudogap. Full microscopic theory of the pseudogap 
state is not our aim here, and in the following we shall deal only with simple modelling of 
appropriate transformation of electronic spectrum and its influence on different physical 
properties, e.g. on superconductivity. 

Considerable simpliflcation of calculations can be achieved if we substitute (jl3|) by 
model interaction of the following form |3H] (similar simpliflcation was flrst used in Ref. 



In fact |T5|) is qualitatively quite similar to ()13|1 and almost do not differ from it quan- 
titatively in most interesting region of |q — Q| < This introduces effective one - 
dimensionality into our problem. 



^In this case all "loop" insertions into interaction lines of perturbation theory and corresponding 
to quantum corrections of higher orders to spin (or charge) fluctuations, and inevitably present in full 
dynamical problem, just vanish. 
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Figure 15: Fermi surface with "hot patches", which are shown by thick hues of the width 
~ C,~^. The angle a determines the size of a "hot patch", a = tt/A corresponds to the 
square Fermi surface. 

Scattering by antiferromagnetic fluctuations in HTSC - oxides is not always most 
intensive at commensurate vector Q = (7r/a,7r/a), in general case Q may correspond to 
incommensurate scattering (see e.g. insert at Fig. |2)). Taking into account the observed 
topology of the Fermi surface with flat parts, as shown in Fig. |2l we can introduce 
another model of scattering by fluctuations of short - range order, which we call the "hot 
patches" model jSU]. In this model we assume the Fermi surface of two - dimensional 
electronic system as shown in Fig. ^1 The size of "hot patches" is determined by 
the angular parameter a. It is well known that flat parts of the Fermi surface usually 
lead to instabilities towards formation of charge (CDW) or spin (SDW) density wave 
and formation of the appropriate long - range order and (dielectric) energy gap at these 
flat parts. We are interested in fluctuation region, when long - range order is not yet 
estabhshed. Fluctuations of the short - range order are again assumed to be static 
and Gaussian, and effective interaction is determined by (|15|). with scattering vectors 
Qx = ^'^.pp, Qy = or Qy = ±2pp, Qx = 0. It is also assumed that fluctuations interact 
only with electrons from these flat ("hot") parts of the Fermi surface, shown in Fig. E[ 
so that this scattering is in fact one - dimensional. In the case of a = 7r/4 we have just a 
square Fermi surface and purely one - dimensional problem. For a < n/A there are also 
"cold" parts of the Fermi surface, where scattering is either absent or small. The choice 
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of scattering vector Q = (±2^^?, 0) or Q = (0, ±2pp) corresponds, in general, to the case 
of incommensurate fluctuations, as the Fermi momentum pp has no direct relation to the 
period of inverse lattice. Commensurate case can also be analyzed within this model [SU] • 
Thus, the main idea of models under discussion reduces to the assumption of strong 
scattering by fluctuations of short - range order, which, according to (fTTj) . (fT3j) . is effective 
in a limited part of momentum space with characteristic size of the order of ^"'^ around 
"hot" spots or patches, which leads to pseudogap transformation of the spectrum in these 
regions. It will be seen in the following, that within our assumptions these models can 
be solved "nearly exactly", and much of the remaining discussion will be devoted to a 
description of this solution. Mostly we shall pay our attention to the discussion of the 
"hot spots" model as more "realistic" and not described in detail in previous reviews. As 
to "hot patches" model - detailed discussion and further references can be found in Ref. 

m 



4 Elementary ("toy") model of the pseudogap. 

Before we go to the analysis of "realistic hot spots model" it may be useful to consider 
an elementary one - dimensional model of the pseudogap, which allows an exact solu- 
tion in analytic form [ST]. Consider an electron moving in one dimension in a random 
field of Gaussian fluctuations with correlation function (in momentum representation and 
identified with an interaction line in appropriate diagram technique) of the following form: 

where k = ^~^{T). The choice of scattering vector Q ~ ±2p^, corresponds to the case of 
incommensurate fluctuations. An exact solution can be obtained in the asymptotic limit 
of ^ ^ CX3 (k — s> 0), i.e. for very large correlation length of fluctuations of short range 
order^. Now we can sum all Feynman diagrams of perturbation theory for an "interaction" 
of the form of ()16|). which in this limit reduces to: 

K//(g) = 27tW^{6{Q - 2pp) + 6{Q + 2pp)} (17) 

Consider the simplest contribution to self - energy of an electron, described by diagram 
shown in Fig. which we write in Matsubara representation: 

dQ 1 2 f°° dx K 1 



ip-^Q J-oo 2tT X^ + lEn + ^p — VpX 



2W^ 



dx K 1 



-oo 2tT {x — iK){x + m) iSn + — VpX 

®Let stress, that this hmit here does not mean the estabhshment of any long - range order. Electron 
moves in the Gaussian random field with special pair correlator, not in periodic system. 
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Figure 16: Simplest diagram for self - energy part of an electron. Wavy line denotes 
interaction Veff{Q). 



i^n + + iVpK 

where, for definiteness, we assumed p ~ +Pf, En > and defined a new integration 
variable x as Q = 2pF + x. We also used here "nesting" property of the spectrum 
ip~2pF = ^^P) which is always vahd in one dimension for the standard = vf{\p\ —Pf)- 
The hmit of ^(T) — > oo (/« ^ 0) should be understood as: 

vpn = vfC^ < Max{27rT, ^p} (19) 



or 

= < 27rT, ^(T) > (20) 

Then (fTHj) reduces to: 

nenP) ^ ——r (21) 

Now, for an "interaction" of the form of ()17j] . there is no problem to write down the 
contribution of an arbitrary diagram of the type shown in Fig. El In such diagram, 
in the n - th order in V^jj(Q) we have 2n vertices, connected, in all possible ways, by 
interaction lines. These lines alternatively^ either "take away" or "bring" the momenta 

^This alternation is important to guarantee that an electron remains close to the Fermi surface (points 
ibpj;^), or large denominators will appear in terms of perturbation theory. This is not important in the case 
of commensurate fluctuations, like period doubling, when we are dealing with tight - binding spectrum 
and "bringing" or "taking away" of any number of momenta Q = (7r/a,7r/a) does not take electron far 
from the Fermi surface. In this case combinatorics of diagrams is different |52j . 
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Figure 17: Diagram of an arbitrary order for the single - electron Green's function. 

Q = 2pp. As a result, the appropriate analytic expression for our diagram contains 
alternating Green's functions \^ (entering n times) and (also entering n times) 

plus an extra (initial) \^ Also we have to take an account of the factor VF^". Finally, 
wc can sec, that contributions of all diagrams in a given order just coincide and their sum 
can be determined from pure combinatorics and is determined just by their number, which 
is equal to n\. It is obvious — there are 2n points (vertices) with "ingoing" or "outgoing" 
interaction lines. Out of these, n points are connected with "outgoing" lines, which can 
"enter" the remaining "free" n vertices in any of n\ ways. Use now the identity 

/"OO rod 1 

j:n\z-^^ dCe-^iCzr^ dCe-^r—- (22) 

n=0 n=0-^° •'^ J- 



Then we can easily sum the whole series for the Green's function and obtain the following 
exact solution: 

= r ""-"'' {iE.r-e^-CW^ =^ G„<^,i,) >c, El = (21 + l)xr (23) 



^°0f course, similar analysis applies to the problem with an arbitrary scattering vextor Q, when we 
have alternating . ^ and -r- — k . We have taken Q = 2pp only to make formula more compact and 

to put the pscudogap precisely at the Fermi level. 
^^In mathematics this is called Borel summation. 
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where we have used the notation: 



(24) 



Now, what has appeared is just the "normal" Green's function of an insulator (of Peierls 
type): 



under the "averaging" procedure of the form: 



< ... > 



d(e~ 



(25) 



(26) 



It is easy to convince yourself (formal proof, as well as many other details on this model 
can be found in Ref. [52]) that ()23|) is just the Green's function of an electron moving in 
an external field of the form 2V cos{2pFX + (p) , with amplitude "fluctuating" according to 
the so called Rayleigh distribution^^: 

2V 

= ^,e-^ (27) 

while the phase (p is distributed homogeneously on the interval from to 27r. 

Performing analytical continuation iei ^ e ± i6 from (j23|) we obtain (for e > 0): 

POO 

•J 

TT 



so that the spectral density 



1 



TC 



{2i 



(29) 



has "non Fermi - liquid like " form, shown in Fig. ^1 Let us stress, that our Green's 

function does not have any poles on the real axis of e, which may correspond to quasi- 
particle energies as required in Fermi - liquid theory. 
Electronic density of states has the following form: 



N{e) 
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dC 



-c 
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exp 



2£_ 

vy2 



Erfi 
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W 



oo 







(30) 



^^This distribution is well known in statistical radiophysics, see e.g.: S.M.Rytov. Introduction to 
Statistical Radiophysics. Part I. "Nauka", Moscow, 1976. 
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Figure 18: Spectral density in the model of pseudogap state: (1) — = 0; (2) — = O.IW; 
(3)— = 0.5W. 
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e/W 



Figure 19: Density of states with pseudogap. 
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Figure 20: Diagram of an arbitrary order for the vertex part of an interaction with an 
external electromagnetic field. 

where Nq[Ef) is the density of states of free electrons at the Fermi level, and Erfi{x) = 
Jq dxe^ is error function of imaginary argument. Characteristic form of this density of 
states is shown in Fig. ^land demonstrate the presence of a "soft" pseudogap around the 
Fermi level. Obviously, this is just the density of states of one - dimensional insulator with 
the energy gap 2V, averaged over fluctuations of this gap with probability distribution 



Remarkable property of this model is the possibility of obtaining an exact solution 
(sum all diagrams) also for the response function to an external electromagnetic field 
[HUEni- The arbitrary diagram for the vertex part, describing the response to an external 
field, can be obtained from an arbitrary diagram for the Green's function (of the type 
shown in Fig. ITTj) by an "insertion" of an external field line into any of electronic lines, 
as shown in Fig. [201 Making such "insertions" into all diagrams of the series fl23|) it is 
possible (after some long, but direct calculations) to sum the whole series for the vertwx 
part and obtain closed expression for response functions, e.g. for polarization operator. 
Details can be found in Refs. [HIl E^l)- However, the structure of an answer is clear 
without any calculations — you must just calculate the response of an insulator with 
fixed gap 2V and afterwards average the result over gap fluctuations with distribution 
function ()27|) . In particular, for polarization operator we obtain the following elegant 
expression (cj^ = 27rmT): 



(EH). 
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+F^w2{enP)F+y2{en + c^mP + q)} =< Il^w^{quj^) >^ (31) 
where automatically appears the product of two "anomalous" Green's functions: 

W* . W 

describing Umklapp processes in a system with long - range order [53j. Due to the absence 
of any long - range order in our model, the value of (jHS)) vanishes after averaging over 
phase, while the average of pair of these functions in two - particle response (j!?T|l is non 
zero. Finally, under the averaging procedure over the gap fluctuations we have simply the 
polarization operator of an insulator (of Peierls type). 

The real part of conductivity for such one - dimensional insulator with fixed gap 2W 
has the following form [H!^ : 



w 



„2 



U > 2W 



Reaw2{id) = \ V w2""^ (33) 
\lj\<2W 

It is seen, that absorption of electromagnetic energy is going through quasiparticle ex- 
citation via the energy gap and is non zero for u > 2W. In the pseudogap state this 
expression must be averaged over fluctuations of W, described by distribution function 
(I2ni) or (EZI). So finally, from §^ we get: 

Rea{u) = 'f^J^dCe-^^^ (34) 

V iw^ ^ 

Appropriate frequency dependence id shown in Fig. |^ We can see characteristic smooth 
maximum of absorption through the pseudogap. 

This elementary model of the pseudogap state is very useful for an analysis of a number 
of problems. It is easily generalized to two - dimensional case for the "hot patches" model 
[50|j . This allows the analysis of the problem of formation of superconducting state on the 
"background" of this (dielectric) pseudogap jSOJ EH ESI EHl EH] • In particular, due to the 
possibility of obtaining an exact solution in closed analytic form, it is possible to study 
rather fine problems of the absence of self - averaging property of superconducting order 
parameter in the random field of pseudogap fluctuations [SI] , showing the possible mech- 
anism of formation of local inhomogeneities ( "superconducting drops" ) at temperatures 
higher than the mean - field critical temperature of superconducting transition. This 
may help to explain e.g. experimentally observed manifestations of superconductivity at 
these high temperatures (like anomalous Nernst effect), which are usually interpreted in 
the spirit of superconducting scenario of pseudogap formation. Possible direct connec- 
tion with the picture of inhomogeneous superconductivity, observed in STM experiments 
[TTIITH], is also obvious. 
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Figure 21: Frequency dependence of the real part of conductivity in the pseudogap state. 
Conductivity is given in units of 4^;^- 

However, advantages of the model determine also its deficiencies. In particular, abso- 
lutely unrealistic is asymptotics of infinite correlation length of pseudogap fluctuations. 
In real systems, as was noted above, this correlation length is usually not larger than few 
interatomic spacings. Besides that, the growth of correlation length inevitably leads to 
the breaking of our assumption of the Gaussian nature of pseudogap fluctuations. Anal- 
ysis of effects of finiteness of correlation length is actually a complicated problem. For 
one - dimensional model such generalization of the model under discussion was proposed 
in Ref. It was shown, that as correlation length ^ becomes smaller, it leads to 

smooth "filling" of the pseudogap, due to the growth of the scattering parameter Vp/C,, 
i.e. of the inverse time of flight of an electron through the region of the size of ~ ^, where 
effectively we have "dielectric" ordering. The method used in Ref. jH^ forms the basis 
of appropriate generalization to two dimensions [IHlllZI; which will be discussed below 
during our analysis of "hot spots" model. Let us also mention the simplified version of 
one - dimensional model with finite correlation length, similar in spirit to the model of 
Ref. [nH EZl 5 proposed in Ref. [HH] and used in Ref. jHSI to analyze the problems of self 
- averaging properties of superconducting order parameter in the "hot patches" model. 



5 "Hot spots" model. 

5.1 "Nearly exact" solution for one — particle Green's function. 

Let us now describe our "nearly exact" solution for "hot spots" model. Consider first 
- order (in Vejf (HHI)) contribution to electron self - energy, corresponding to simplest 
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Figure 22: Diagrams of second order in effective interaction with pseudogap fluctuations, 
diagram, shown in Fig. ^1 

S(5„,p) = Y: K//(q) . _\ (35) 

For large enough correlation lengths ^, the main contribution to the sum over q comes 
from the vicinity of Q = (vr/a, vr/a). Then we can write: 

■Cp+q = 'Cp+Q+k ~ 'Cp+Q + Vp+qk (36) 

where Vp+q = is the appropriate velocity of a quasiparticle on the Fermi surface. 

Then fl35|) is easily calculated and we get: 

^i^nP) = -. T — r-y^ ^ : (37) 

ISn - ^p+Q + HI^p+qI + \Vp+Q\) f^StgUEn 

where k = Let us stress that both here and below "linearization" of the quasiparticle 
spectrum (jHUI) under the integral (jH^j) is performed only over small, due to large enough 
^, correction (of the order of vfC,~^) to the spectrum close to the Fermi surface, while the 
forms of the spectrum itself and ^p+q are given by the general expression (fTUI) with 

Consider now correction of the second order, shown in Fig. 1221 Using H15|l we obtain: 
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2 I h'^ t^^ 4- y'^ _l_ h'^ i<-2 I ^.2 



iSn - ^p+Q - t'p+qfcla; " Wp+Q^iy ie„ - Cp - ^'p(fclx + ^23:) " ^^p(fcly + ^2,j) 

1 

i^n — ^p+Q — ^p+Q^la; ^ ^p+Q^ly 
dki /■ (ik2 K K K K 



(38) 



7r2 7 772^2 + K? + fc? ^2 + fcf^ ^2 + fc| 



1 1 
- Cp+Q - ^^p+Q^la; - ^p+Q^iy ^ ^ t^p(fclx + k2x) - Vpihy + k2y) 

1 

- ^p+Q - 'fp+Qfc2a; - I'p+Qfc2y 



(39) 



where we have used the exphcit form of the spectrum H10() . from which it follows, in particular, that 
^p+2Q = £.p, Vp+2Q = Vp for Q = (7r/a,7r/a). If the signs of Vp and Vp+Q, as weh as of and Wp+q 
coincide, the integrals in 1)38(1 and ((39|l are completely determined by contributions from the poles of 
Lorentzians, which describe interaction with fluctuations of short - range order, so that after elementary 
contour integration we get: 

^^"^ " ^^^^ " [^en Cp+Q + ^{\vUci\ + K+Q\h? - Cp + *2(|^;-| + Iv^l)^ ^^^^ 

Here and in the following, for definiteness, we assume e„ > 0. 

It is not difficult to convince yourself, that in case of coinciding signs of velocity projections at "hot 
spots" , similar calculation is valid for an arbitrary diagram of higher order. 

Thus, in case of coinciding signs of velocity projections at the Fermi surface Vp and 

"fp+Q, and those of and "Wp+q, Feynman integrals in any diagram of arbitrary order 

are determined only by contributions form the poles of Lorentzians in ()15j) and are easily 

calculated. Similar situation holds also in the case of velocities at "hot spots" , connected 

by vector Q, are perpendicular to each other. In this case the contribution of an arbitrary 

diagram of - th order in (|T3j) for electron self - energy has the following form: 

2N-1 -1 

SW(.„p) = W'' n -■ TTT-- (41) 

where = ^p+q and vj = \v^+q\ + \vI+q\ for odd j and ^^(p) = ^p and Vj = \v^\ + \vy\ 
for even j. Here nj is the number of interaction lines, surrounding j - th Green's function 
(counting from the first one in diagram) and we again take e„ > 0. 

In Ref. jlH] we have studied in detail when these conditions on velocity projections are 
satisfied in the points of the Fermi surface connected by vector Q ( "hot spots" ) and have 
presented explicit examples of appropriate geometries of the fermi surface, which can be 
realized for specific relations between parameters t and t' in (jiup . In these cases expression 
()41|) is exact, with only limitation being related to our use of "linearization" (jHUj) . In all 
other cases (for other relations between t and t') we use fl4ip as rather successful Ansatz 
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for the contribution of an arbitrary order, obtained by simple continuation of spectrum 
parameters t and t' to the region of interest to us. Even in most inappropriate one - 
dimensional case ^7\, corresponding to square Fermi surface, appearing for (fTUIl if t' = 
and fi = 0, the use of this Ansatz produces results (e.g. for the density of state) which 
are very close quantitatively jHOl to the results of an exact numerical simulation of this 
problem |tilj . 

Generally speaking, note that in standard diagrammatic approaches we usually per- 
form summation some subseries of diagrams, which are considered as dominant over some 
smallness parameter. Here we are dealing with much more rare situation — we can sum 
the whole diagram series, though contribution of each diagram is calculated probably 
approximately. It is in this sense that we use the term "nearly exact" solution. 

Using Ansatz (|4ip we can see, that the contribution of an arbitrary diagram with in- 
tersecting interaction lines is actually equal to the contribution of some diagram of the 
same order without intersections of these lines ^7\. Thus, in fact we can limit ourselves to 
consideration of only diagrams without intersecting interaction lines, taking the contribu- 
tion of diagrams with intersections into account with the help of additional combinatorial 
factors, which are attributed to "initial" vertices or just interaction lines jHZI- As a result 
we obtain the following recursion relation (continuous fraction representation ^7\), which 
gives an effective algorithm for numerical computations 

1 



i^n - ^kip) + ikvkK - Sfe+i(£:„^p) 

= {Gofe(^nep)-Sfc+l(5„^p)}"' (42) 

s.(^nep) = W', .y ^ T ( ,^ (43) 

- a(p) + ^kvk^^ - Efc+i(£:„4p) 
Graphically this recursion relation for the Green's function is shown in Fig. 1221 The 
"physical" Green's function is obtained as: G{eni-p) = Gk=o{snC,p)- In we also defined: 

Gokisn^p) = -. , (44) 

^^n - 4fc(P) +lkVkK 

Combinatorial factor: 

v{k) = k (45) 

for the case of commensurate fluctuations with Q = (7r/a,7r/a) ^E7\, if we don not take 
into spin structure od interaction (CDW - type fluctuations). For incommensurate CDW 
fluctuations [IFfj : 

,(k) = ( f * (46) 

^ ' [ I for even k ^ ' 

If we take into account (Heisenberg) spin structure of interaction with pseudogap 
fluctuations in "nearly antiferromagnetic Fermi - liquid" (spin - fermion model [47 ), 
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Figure 23: Graphical representation of recursion relation for the single - particle Green's 
function. 

combinatorics of diagrams becomes more complicated. Thus, spin - conserving scattering 
processes obey commensurate combinatorics, while spin - flip scattering is described by 
diagrams of incommensurate case ( "charged" random field in terms of Ref . [47] ) . In this 
model recursion relation for the Green's function again is given by (j43|) . but combinatorial 
factor v{k) takes now the following form |47j : 

(U) - j ^ '^^^ ^ (A7\ 

1 g for even k 

The obtained solution for the single - particle Green's function is asymptotically exact 
in the limit of ^ — > oo, when solution can be found also in analytical form 07]. It 
is also exact in trivial limit of ^ — 0, when for fixed values of W interaction (jl5|) just 
vanishes. For all intermediate values of ^ our solution gives, as already noted, very 
good interpolation, being practically exact for certain geometries of the Fermi surface, 
appearing for certain relations between parameters of the spectrum (jlUp [15] . Note also 
that our formalism can be easily used also to describe pseudogap within superconducting 
scenario ^5] jHSj- Our preference for "dielectric" scenario is based mainly on physical 
considerations. 

Using ()42j) we can easily perform numerical calculations of single - electron spectral 
density: 

A{Ep) = --ImG^iEp) (48) 
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Figure 24: Spectral density in the "hot spots" model, the case of incommensurate fluctua- 
tions, Ka = 0.01: (1) - at the "hot spot" Pxa/n = 0.142, p^a/vr = 0.587, (2) - close to "hot 
spot" atp^a/iT = 0.145, Pyo/vr = 0.843, (3) -far from "hot spot" aipxa/n = Py/71 = 0.375. 

which can also be determined from ARPES experiments IJ^. In fl48j) G^{Ep) represents 
retarded Green's function obtained by the usual analytic continuation of from Mat- 
subara frequencies to the real axis of E. Analogously we cam compute single - particle 
density of states as: 

N{E) = Y: AiEp) = -- E IrnG'^iEp) (49) 
p ^ p 

Details of these calculations and discussion of the results for our two - dimensional model 
can be found in Refs. ^ZlEHj- Here we shall demonstrate only few most important results. 

As a typical example in Fig. |^ we show the results for spectral density of 
electrons for the incommensurate (CDW) case. We can see that spectral density close to 
the "hot spot" has an expected non Fermi - liquid like form and there are no well defined 
quasiparticle peaks. Far from the "hot spot" spectral density is characterized by a narrow 
peak, corresponding to well defined quasiparticles (Fermi - liquid). In Fig. EHlfrom Ref. 
|T7j we show the product of the Fermi distribution and spectral density at different points 
of the "renormalized" Fermi surface, defined by the equation Sp — ReTj{E = Op) = 0, 
where the "bare" spectrum is given by (fTUj) with t = — 0.25eV,t' = — 0.35t and for 
hole concentration Uh = 0.16, with coupling constant in (fTT|) g = 0.8eV^ and correlation 
length ^ = 3a (commensurate case, spin - fermion model). We can clearly see complete 
qualitative agreement with ARPES data discussed above with quite different behavior 
close and far from the "hot spot" . Finally, in Fig. 1211 taken from Ref. we show 
calculated (for spin fermion - model) positions of the maximum of A{ujk.) for two different 
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Figure 25: The product of spectral density and Fermi distribution function at different 
points on the Fermi surface, shown at the insert. Spin - fermion model, correlation length 
^ = 3a. 




Figure 26: Positions of the maximum of the spectral density for different values of correla- 
tion length ^ and hole concentrations, calculated for spin - fermion model and compared 
with ARPES data for Bi2Sr2Cai-xDyxCu20s+& with x = 1 (triangles) and x = 0.175 
(diamonds). 
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Figure 27: Single - electron density of states (in units of l/ta?) for for different combina- 
torics of diagrams and t' /t = —0.4, fi/t = —1.3, which is typical for HTSC - cuprates: (1) 
— incommensurate case, (2) — commensurate case, (3) — combinatorics of spin - fermion 
model, (4) — in the absence of pdeudogap fluctuations. We assume W/t = 1 and corre- 
lation length na = 0.1. At the insert: Density of states for commensurate combinatorics 
and: (1)— /ta = 0.1; (2)— /ta = 0.01 

concentrations of holes compared with appropriate experimental ARPES data of Ref. [^21 
for Bi2Sr2Cai_xDyxCu20s+s- The point is, that positions of the maxima of spectral 
density in the plane of (co", k), determined from ARPES, in an ideal system of the Fermi - 
liquid type define dispersion (spectrum) of quasiparticles (cf. Fig. Efa)). For overdoped 
system the values of = 0.22 and ^ = a were assumed during these calculations. The 
results obtained demonstrate rather well defined dispersion curves both along the diagonal 
of Brillouin zone, and in the direction (0,0) — (vr, 0). For underdoped system it was 
assumed that rih = 0.16 and ^ = 3a. In this case, in diagonal direction we again see 
spectral curve crossing the Fermi level, while close to "hot spots" (in the vicinity of 
(tt, 0)) there is seen only a smeared maximum of spectral density remaining approximately 
200 meV below the Fermi level (pseudogap). In general, agreement between theory and 
experiment is rather satisfactory. 

Ley us now consider the single - electron density of states, defined by the integral of 
the spectral density A{Ep) over the whole Brillouin zone. Detailed calculations of the 
density of states in "hot spots" model were performed in Ref. |1H]- As an example, in 
Fig. 123 we show appropriate data for the Fermi surface topology typical for HTSC - 
systems. We can see, that for typical value t' /t = —0.4 there is a shallow minimum in the 
density of states (pseudogap), which is only slightly dependent on the value of correlation 
length ^. At the same time, e.g. for t' /t = —0.6 (which is untypical for HTSC - cuprates) 
there are "hot spots" on the Fermi surface, but pseudogap in the density of states is 
practically unobservable |1H] • We can see only smearing of Van Hove singularity, which is 
present for an ideal case in the absence of pseudogap scattering. In this sense most clear 
manifestations of the pseudogap behavior are not in the density of states, but in spectral 



43 



v(i) «(2) m 



qra 



Figure 28: General form of higher order correction to the vertex part, 
density, which is in general accordance with experiments. 

5.2 Recurrence relations for the vertex part and optical con- 
ductivity. 

To calculate optical conductivity we need to know the vertex part, describing electromag- 
netic response of our system. This vertex part can be found using the method, proposed 
for similar one - dimensional model in Refs. [^3 IHS] • Below we are following Ref . . 
Above we have seen, that any diagram for irreducible vertex part can be obtained by the 
insertion of an external field line in appropriate diagram for electron self - energy [HI] . 
As in our model it is sufficient to take account only of self - energy diagrams without in- 
tersecting interaction lines with additional combinatorial factors v{k) in "initial" vertices, 
to calculate vertex corrections it is possible to limit ourselves only to diagrams of the 
type shown in Fig. |2H1 Diagrams with intersecting interaction lines are also accounted 
for automatically. Now we immediately obtain the system of recursion equations for the 
vertex parts, shown graphically in Fig. EHl To obtain appropriate analytic expressions, 
consider the simplest vertex correction, shown in Fig. EDI (a). Performing calculations for 
T = in RA - channel, we can easily obtain its contribution as: 



r(l)RA 

K 

W^G^oie,^i{p)+iviK)-G^oie + Lo,^i{p + q)-iviK)} 



1 



^ + 6(p)-a(p + q) 
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Figure 29: Recursion relations for the vertex part. 
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Figure 30: Simplest corrections to the vertex parts. 
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^ + ei(p) -6(p + q) J 

- «'%.V^,«P)G?.(. +^.«..,) {l + (50) 

where, during calculations of integrals, we have used the following identity, valid for free 
- particle Green's functions: 

GiieQG^,iE + ^^p+q) = {G^oH^) - 0^,(6 + ^^p+j} T^—rr (51) 



"Dressing" internal electronic lines we obtain diagram, shown in Fig. 120 (b), and using 
the identity: 



G^{e^^)G^{e + cu^p+q) = {G^H^) - G^'ie + o^^p+q)} 

1 



^ - ^P+q + ^p - Sf (£ + ^^p+q) + ^iHp 

valid for full Green's functions, we can write down the contribution of this diagram as: 
J«^(ep; e + Lup + q) = WMl)Gti£, Cp)G?ie + Cp+q) |l + 

2iVlK \ 



(52) 



ii (P + q) + ii (P) - (e + c^Cp+q) + (eCp) 



Jf^(ep;e + tjp + q) 

(53) 



Here we assumed, that interaction line in the vertex correction diagram of Fig. IHO (b) 
"transforms" self - energies Sf'"^ of internal electronic lines into S^'"^, in accordance with 

our main approximation for self - energies used above (cf. Fig. I^Bj)^^. 

Now it is not difficult to write down similar expression for the diagram of general form, 
shown in Fig. 1201 (c): 

J^^iep- £ + cop + q) = W\{k)Gt{e, (p)G^ie + cj, ^p+q) 1 1 + 

+ - 6(P + q) + a(p) - + c^fp+q) + ^tMp) ] ' ^^P'" + "P + ^^ 

(54) 

Accordingly, we can write dowm the following fundamental recursion relation for the 
vertex part, shown in Fig. 1221 

J^A (ep; e + cop + q) - 1 + W\{k)Gt{e, Cp)G^{e + co, ^p+q) | 1 + 



^■^This assumption is justified by the fact, that it guarantees vahdity of certain exact relation, following 
from the Ward identity (see below) |64j . 
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Figure 31: Diagrammatic representation for two - particle function ^^^{qu). 

2ivknk 

+ ^ - 6(P + q) + efe(p) - + ^^P+q) + ^t+Mp) 

(55) 

"Physical" vertex J^^{ep; e + up + q) is defined as J^(ep; e + tup + q). Recursion 
procedure ()55p takes into account all diagrams of perturbation theory for the vertex part. 
For K — > — > oo) (fH^ redoces to the series, studied in Refs. jHI] (see also Ref. 
|T7j). which can be summed exactly in analytic form. Standard "ladder" approximation 
is obtained from our scheme putting all combinatorial factors v{k) in (j55|) to unity [HH] . 
Conductivity is expressed via retarded density - density response function x^(gu;) [U7] : 

a{uj) = eHim[-^]x''{qio) (56) 

where e is electronic charge, 

X'^iquj) = CO {$^^(Ogcu) - <l>^^(00cu)} (57) 
where two - particle Green's function ^^^{equj) is defined by loop diagram shown in Fig. 

EU 

Direct numerical computations confirm, that the recurrence procedure (j55p satisfies 
an exact relation, which directly follows (for — > 0) from the Ward identity ^7\: 

$^^(00a;) = -^^^ (58) 

UJ 
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J;f^(ep;e + wp + q) 



where N{Ep) is the density of states at the Fermi level Ep = /i. Actually, this is our 
main motivation for the Ansatz, used during the derivation of (j53|) . (j54p and (j55|) . 

Finally, conductivity is written in the following symmetrized form, convenient for 
numerical computations: 




(59) 

where we also accounted for an extra factor of 2, due to summation over spin. 

Direct numerical calculations |^ were performed, for different values of parameters 
of the "bare" spectrum (|Tn|l . using ([Hlljl . with recursion procedure starting at 

some high value of k, where all and Jk were assumed to be zero. Integration in (jSHI) was 
made over the Brillouin zone. Integration momenta are naturally made dimensionless with 
the help of lattice constant a, all energies below are given in units of transfer integral t. 
Conductivity is measured in units of universal conductivity of two - dimensional system: 
(To = ^ = 2.5 10"^ Ohm~^, and density of states — in units of 1/ta^. For definiteness, 
below we always take W = t. 

First let us consider Fermi surfaces close to the case of half - filled band fi = and t' = 
0, shown (in the first quadrant of Brillouin zone) in Fig. EH (a). We know that for /i = 
and t' = Fermi surface is simply square (complete "nesting"), so that we practically 
are dealing with one - dimensional case, analyzed long ago in Refs. jSH EH IHSl- Results 
of our calculations for real part of optical conductivity in our two - dimensional model, 
for the case of spin - fermion diagram combinatorics and different values of correlation 
length of AFM short - range order (parameter k = where ^ is measured in units 
of lattice constant a) are shown in Fig. ESI Qualitative behavior of conductivity is 
quite similar to those found for one - dimensional model (for the case of incommensurate 
fluctuations of CDW - type) in Refs. [HlllSni- It is characterized by the presence of large 
pseudogap absorption maximum (appropriate densities of states with pseudogap close to 
the Fermi level are shown at the insert in Fig. IH^ at a; ~ 2W and also by a maximum 
at small frequencies, connected with carrier localization in static (in our approximation) 
random field of AFM fluctuations. Localization nature of this maximum is confirmed by 
its transformation into characteristic "Drude - like" peak (with maximum at = 0) if 
we perform calculations in "ladder" approximation, when combinatorial factors v{k) = 1, 
which corresponds to "switching off" the contribution from diagrams with intersecting 
interaction lines, leading to two - dimensional Anderson localization |^ l68j. Qualitative 
form of conductivity in this case is also quite similar to those found in Ref . [HH] • Narrowing 
of localization peak with diminishing correlation length of fluctuations can be explained, 
as was noted in Ref. jSH]; by suppression of effective interaction (|15|) for small ^ (with 
fixed value of W), leading to general suppression of scattering, also on the "cold" part 
of the Fermi surface. Note that general behavior of the density of states and optical 
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Figure 32: Fermi surfaces for different values of t' and chemical potential ji. (a) corre- 
sponds to = and the following values of t' /t: - I; -0.2 - II; -0.4 - III; -0.6 - IV. (b) 

corresponds to t' = and the following values of /x/t: - I; -0.3 - II; -0.5 - III; -0.6 - IV. 
At the inserts we show energy dependences of spectral densities for spin - fermion model 
for fi;a = 0.1 at points of the momentum space denotes by stars. 
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Figure 33: Real part of optical conductivity in spin - fermion model for the case of square 
Fermi surface [fi = 0, t' = 0) for different values of inverse correlation length Ka: 0.1 - 1; 
0.2 - 2; 0.5 - 3. At the insert we show appropriate densities of states. 

conductivity obtained here is in complete qualitative agreement with the results, obtained 
for the similar two - dimensional model of Peierls transition by quantum Monte - Carlo 
calculations in Ref. 

Now, if we keep /i = and "switch on" the transfer integral t' between second nearest 
neighbors in ()10|) . we shall obtain Fermi surfaces different from square one, as shown in 
Fig. 1221 (a)- At the insert on this figure we also show the energy dependence of spectral 
density PH|) at several typical points on these Fermi surfaces. It is seen that spectral 
density demonstrates characteristic "non Fermi - liquid" behavior, of the type studied in 
Refs. I3H] , practically everywhere on the Fermi surface, until this surface is not very 
different from the square, despite the fact that the "hot spot" in this case lies precisely 
at the intersection of the Fermi surface with diagonal of the Brillouin zone. Appropriate 
dependences of the real part of optical conductivity are shown in Fig. |^ At the insert 
in this figure we show appropriate densities of states. It is seen that as we move away 
from the complete "nesting", pseudogap absorption maximum becomes more shallow, 
while localization peak (in accordance with the general sum rule for conductivity) grows. 
Note, however, that pseudogap absorption remains noticeable even when pseudogap in 
the density of states is practically invisible (curves 4 in Fig. IMjl . 

Let us return now to the case of t' = and change the value of /i, so that Fermi 
surfaces are close to the square, as shown in Fig. 1221(b). Strictly speaking, "hot spots" 
on these surfaces are absent, but spectral density shown at the insert on Fig. 1221(b), is still 
typically pseudogap like. Appropriate dependences of the real part of optical conductivity 
are shown in Fig. (23 
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Figure 34: Real part of optical conductivity in spin - fermion model for jj, = and 
Ka = 0.1 and for different Fermi surfaces, obtained from the square after "switching on" 
the transfer integral t'/t: - 1; -0.2 - 2; -0.4 - 3; -0.6 - 4. At the insert — appropriate 
densities of states. 
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Figure 35: Real part of optical conductivity in spin - fermion model for t' = and 
= 0.1 and for different Fermi surfaces, obtained from square as we move from the case 
of half - filled band. Chemical potential /i/t: 0-1; -0.3 - 2; -0.5 ~ 3; -0.6 - 4. 
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Figure 36: Real part of optical conductivity for t' /t = —0.4 and fi = for Ka = 0.1 and for 
different combinatorics of diagrams: 1 - spin - fermion combinatorics; 2 - commensurate 
case. Dashed curve — "ladder" approximation. At the insert — appropriate densities of 
states. 

Consider now geometry of the Fermi surface with "hor spots" typical for most of HTSC 
- oxides, like that shown in Fig. EI First, in Fig. [HHlwe show the real part of optical 
conductivity, calculated (for different combinatorics of diagrams) for characteristic value 
of t' = —OAt and for chemical potential /i = 0, when "hot spots" are on the diagonals 
of the Brillouin zone. It is seen that pseudogap behavior of conductivity is conserved 
even in the case of practically absent pseudogap in the density of states (shown at the 
insert in Fig. Dashed curve in Fig. IHHl shows the result of "ladder" approximation, 
demonstrating disappearance of two - dimensional localization. In general, as correlation 
length of the short - range order diminishes, we observe "smearing" of the pseudogap 
maximum in conductivity. 

For most copper oxide superconductors characteristic geometry of the Fermi surface 
can be described by t' = —OAt and fx = —1.3t fT?. In this case, results of our calculations 
of optical conductivity for different values of inverse correlation length k are shown in 
Fig. EH (for the case of spin - fermion combinatorics), here we introduced additional 
weak scattering due to inelastic processes via standard replacement u ^ u + i'j [70] , 
which leads to the appearance of a narrow "Drude - like" peak for < 7 (destruction 
of two - dimensional localization due to phase decoherence) . It is easy to see that with 
the growth of inelastic scattering rate 7, localization peak is completely "smeared" and 
transformed to the "usual" Drude - like peak at small frequencies. Pseudogap absorption 
maximum becomes more pronounced with the growth of correlation length ^ (diminishing 
/t). Similar results for the "hot patches" model were obtained in Ref. [7T] 
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Figure 37: Real part of optical conductivity in spin - fermion model for t' /t = —0,4 and 
fi/t = —1.3 and for different values of correlation length Ka: 0.05 - 1; 0.1 - 2; 0.2 - 3. 
Damping factor j/t = 0.005. 

Direct correspondence of these results with experimental data shown in Fig. |3] PH] is 
obvious. Simple estimates show that characteristic values of conductivity in theory and 
experiment are also of the same order of magnitude. In principle, it is quite possible to 
make quantitative fit to the experiment, varying parameters of our model. 

5.3 Interaction vertex for superconducting fluctuations. 

Now we are going to discuss superconductivity formation on the "background" of pseu- 
dogap fiuctuations. To take into account pseudogap fiuctuations during the analysis of 
Cooper instability and derivation of Ginzburg - Landau expansion we need knowledge of 
the vertex part, describing electron interaction with an arbitrary fiuctuation of supercon- 
ducting order parameter (gap) of a given symmetry [72j: 

A(p,q) = Aqe(p) (60) 

where the symmetry factor determining the type (symmetry) of pairing is: 

^ r 1 ( s - wave pairing) 

^ ' 1 cos Ps^a — cos Pytt ( d^2_y2 - pairing) 

and we always assume singlet pairing. 

Almost immediately we can write down recursion relations [721 ES] for "triangular" 
vertices in Cooper channel, similar to those introduced above during the analysis of re- 
sponse to electromagnetic field |66j- The vertex part of interest to us can be defined 
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as: 



Tier 



-p + q) 



-e„,qje(pj 



(62) 



Then rp(e„, — £„, q) is determined by the recursion procedure of the following form: 



Tpk-i{en, ~en,q) = 1 ± W"V(fc)G'fc(e„, p + q)G'fc(-£„, p)|l + 

2ikK,Vk 



G'fc^(£„,p + q) - Gfc^ (-£„,, p) - 2zA;Kt;fe 



>rpA:(e^ 



(63) 



which is represented by diagrams shown in Fig. |381 "Physical" vertex corresponds to 
Tpk=o{£n, —£n,<i)- Additional combinatorial factor r{k) = v{k) for the simplest case 
of charge (or Ising like spin) pseudogap fluctuations analyzed in Ref. [Z2j. For the 
most interesting case of Heisenberg spin (SDW) fluctuations, which are mainly considered 
below, this factor is given by jUl US] : 



r{k) 



k for even k 
^ for odd A; 



(64) 



The choice of the sign before W"^ in the r. h. s. of (jU^ depends on the symmetry of su- 
perconducting order parameter and the type of pseudogap fluctuations [Z2lEn]- Summary 
of all variants is given in Table I. 

Table I. The choice of sign in recursion procedure for the vertex part. 



Pairing 


CDW 


SDW (Ising) 


SDW (Heisenberg) 


s 


+ 




+ 


d 




+ 





From this table we can see, that in most interesting case of d - wave pairing and Heisen- 
berg pseudogap fluctuations this sign is "— ", so that we have recursion procedure with 
alternating signs. At the same time, for the case of s - wave pairing and the same type of 
fluctuations we have to take this sign and signs in recursion procedure are always the 
same. In Ref. [72] it was shown that this difference in types of recursion procedure leads 
to two different variants of qualitatively different behavior of all the main characteristics 
of superconductors. 



5.4 Impurity scattering. 

Scattering by normal (nonmagnetic) impurities is easily taken into account in self - con- 
sistent Born approximation, writing down "Dyson's equation" shown diagrammatically in 
Fig. inni(a), where in addition to Fig. EHl we have added impurity scattering contribution 
to electron self - energy. 
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Figure 38: Recursion equations for "triangular" vertex in Cooper channel. 
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Figure 39: Recursion equations for the Green's function (a) and "triangular" vertex (b) 
with the account of scattering by impurities. 
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As a result, recursion equation for the Green's function is written as: 



^'^^"P^ Got (^nP) - pU^ Ep G(e.p) - WH{k + l)Gk+i{enP) ^^^^ 

where p is the concentration of impurities with point - hke potential U , and impurity 
scattering self - energy contains full Green's function G(£„p) = Gk=o{SnV)^ which is to 
be determined self - consistently by our recursion procedure. A contribution to impu- 
rity scattering self - energy form the real part of Green's function reduces, as usual, to 
irrelevant renormalization of the chemical potential, so that (j65|) reduces to: 

^'^^"""^ ^ i{en - pU' Ep ImGisnP) + kv^n) - Uv) - W^v{k + l)Gfc+i(e„p) ^^^^ 
Thus, in comparison with impurity free case, we have just a substitution (renormalization): 

£n £n - pW^ ImG{enP) = SnVe (67) 
P 

of/2 

Ve = l-^T.^^G{enP) (68) 

p 

If we do not perform fully self - consistent calculations of impurity self - energy, in the 
simplest approximation we simply have: 

Sn- pU^ ImGooienP) = SnVe = £n + josignEn (69) 
p 

= 1 + ri (^0) 

I 

where 70 = npU'^No{0) is the standard Born impurity scattering rate (Ai'o(O) is the density 
of states of free electrons at the Fermi level). 

For "triangular" vertices of interest to us, recursion equations with the account of 
impurity scattering is shown diagrammatically in Fig. [HHl (b) . For the vertex, describing 
electron interacting with superconducting order parameter fluctuation ()60p with d - wave 
symmetry ()61|). this equation is considerably simplified, as the contribution of the second 
diagram in the r. h. s. of Fig. [HHl (b) in fact vanishes due to J2p e(p) = (cf. discussion of 
similar situation in Ref. jSOl)- Then the recursion equation for the vertex takes the form 
(lUH)) . where Gfc(±£„p) are given by (jHSI), i-e. are just Green's functions "dressed" 
by impurity scattering, determined gy diagrams of Fig. EHl (a)- For the vertex describ- 
ing interaction with order parameter fluctuations with s - wave symmetry, we have the 
following equation: 

rpfc-i(£:„, -Sn, q) = 1 + pf/^ Y '^(^«' P + ^)Gi-£n, P)rp(e, -€n, q) ± 



58 



±W\{k)Gk{en, p + q)Gk{-en, p) U + 



(71) 



p + q) - p) - 2ikKVk 



2ikKVk 




where Gfc(±e„p) is again given by §B>, dSEI), and the sign before is determined 
according to the rules formulated above. The difference with the case of the vertex, 
describing interaction with d - wave fluctuations, is in the appearance of the second term 
in the r. h. s. of (f7T|) . so that we have to substitute: 



Now the self - consistent procedure look as follows. We start from "zeroth" approx- 
imation G = Gqo, Fp = 1, then in Eqs. (jUBj) . (fTTj) we just have 77^ = r^r = 1 — 
pW^/enJ^pI^Goo^EnP)- Then perform recursions (starting from some big enough value 
of k) and determine new values of G = Gk=Q and Fp = F^^o- Again calculate 77^ and ?7r 
using and (f72|) . put these into (jHSI), (f7T|) etc., until convergence. 

While analyzing vertices with d - wave symmetry, we have simply to put rjY = 1 aX all 
stages of calculation. In fact, in these case there is no serious need to perform fully self 
- consistent calculation, because it leads only to insignificant corrections to the results of 
non self - consistent calculation, using only simplest substitution (jHEI) [13 • 

As an illustration, in Fig. 1301 we show comparison of ARPES data of Ref. for mo- 
mentum dependence of a = ImIl{E = 0, p), taken from Fig. ^l(c), with the results of non 
self - consistent calculation using (jHH), (jZOl)- Assumed values of parameters, typical 
for HTSC, are shown on the figure, for the parameters of the spectrum ()10|) we have taken 
t = 0.25eV, t' = —OAt, while chemical potential was calculated for two limiting doping 
levels. It can be seen that we obtain correct order of magnitude estimate of anisotropy 
in momentum space, but the general form of this dependence is only qualitatively similar 
to that observed in the experiment, which lies in between two calculated curves. More 
or less similar results are obtained also in the case of spin - fermion combinatorics. In 
principle, rather approximate agreement of calculated results with experiment is not very 
surprising. Our model is certainly oversimplified, and experimental data are also not 
very precise. Besides, we practically know nothing about the values of parameters of the 
model, appropriate for the system, studied in these experiments. 

5.5 Superconducting transition temperature and Ginzburg — 
Landau expansion. 

Critical temperature of superconducting transition is determined by the equation for 
Cooper instability of the normal phase: 



1 ^ ?7r = 1 + pU^ G{£n, P + q)G{-en, p)Fp(£:, q) 



(72) 



p 



l-yx(0;T)=0 



(73) 
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Figure 40: Comparison of ARPES data for the static contribution to the imaginary part 
of electron self - energy a = ImT,{E = 0, p) with calculated results for the "hot spots" 
model with impurity scattering, for typical set of model parameters, shown in the figure. 
Shown are the results for commensurate combinatorics and two values of doping: 0% and 
20%. Momentum is given in units of ^, and its change between and 1 corresponds 
to measurement point moving along the Fermi surface from the diagonal of the Brillouin 
zone towards the vicinity of the point (0, tt). 



where generalized Cooper susceptibility is defined by diagram shown in Fig. |^ and is 
equal to: 

x(q;T) = -T^^G(£„p + q)G(-£„,-p)e2(p)rp(£„,-£„,q) (74) 

s„ p 

Pairing coupling constant V is assumed to be nonzero in some layer of the width of 
2uc around the Fermi level and determines the "bare" transition temperature Tco in the 
absence of the pseudogap fluctuations via standard BCS equation^^: 



dpx / dpy—- — - (75) 



2VT ^ f^h , r^h , e2(p 

where m = [3^^] is dimensionless cut - off parameter of the sum over Matsubara frequen- 
cies. All calculations [721 ESI were performed for the typical spectrum of quasiparticles 
(fTUI) . for different relations between t, t' and /i. Choosing, rather arbitrarily, ujc = OAt 
and Tco = O.Olt, we can easily find the appropriate value of pairing interaction V in 
fl75|l . leading to this given value of Tco for different types of pairing. In particular, for 
t' /t = —0.4, jJi/t = —1.3 and for s - wave pairing we get ^ = 1, while for dx^^-y^ - pairing 
we obtain = 0.55. 

To determine Tc we need only the knowledge of Cooper susceptibility for q = which 
considerably simplifies all calculations In general case, e.g. to derive coefficients of 
Ginzburg - Landau expansion we need to know x(g; T) for arbitrary (small) q. 

^"^We do not discuss the microscopic nature of this interaction - it can be due to an exchange by AFM 
spin fluctuations, phonons, or combination of both. 



60 



En P*1 




Figure 41: Diagrammatic representation for generalized susceptibility x(q) in Cooper 
channel. 

Ginzburg - Landau expansion for the difference of free energies of superconducting 
and normal states is written in the following standard form: 

- F„ = A|Aq|2 + Aq|2 + ||Aq|^ (76) 

and is determined by the loop expansion of free energy of an electron in the field of 
fluctuations of suoerconducting order parameter ()60|) . shown in Fig. 021 

It is convenient to normalize Ginzburg - Landau coefficients A, B, C by their values 
in the absence of pseudogap fluctuations, writing the following representation [^2] : 



A = AoKa\ C = CoKc; B = BoKb, (77) 

^0 = ^o(O)^^ < e'ip) >; Co = A^o(O)^^ < |v(p)|V(p) >; 

Bo = ^o(O)^ < e\p) >, (78) 

where the angular brackets denote the usual averaging over the Fermi surface: < . . . >= 
T^Q) Sp'^(Cp) • • •) where A^'o(O) - density of states on the Fermi level for free electrons. 
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Then we obtain the following general expressions [7^ : 



x(0;T)-x(0;T,) 
J^A = z (79) 

^0 

irc = lim ^^q''^-^-^^°''^-^ (80) 
Kb = ^Y.T. e\v){G{enV)G{-en, -p))'(rp(e„, -e., 0))^ (81) 

^0 p 



which were used for direct numerical calculations. In the presence of impurities all Green's 
functions and vertices, entering these expressions, should be calculated using Eqs. 
and (f7T|). derived above. 

The knowledge of the coefficients of Ginzburg - Landau expansion allows to determi- 
nation of all major characteristics of a superconductor close to the transition temperature 
Tc. Coherence length is defined as: 

em _Ko 



eBcs(T) 

where ^bcs(T) is the value of this length in the absence of the pseudogap. Penetration 
depth is: 

where we again normalized to the value of XscsiT) in the absence of pseudogap fluctua- 
tions. Analogously, normalized slope of the upper critical field close to is: 



dHc2 




dT 




dHr.2 




dT 





J-cO i^c 



Specific heat discontinuity at the transition: 



{Cs — Cn)TcO ^cO Kb 

Results of these calculations for the cases of charge (CDW) and Ising like spin (SDW) 
fluctuations of short range order can be found in Ref. |72j- Here we shall concentrate 
mainly on the analysis of most important and interesting case of Heisenberg spin (SDW) 
fluctuations and also on the discussion of the role of impurity scattering (disordering) 
[7^ [73j. Due to particular importance of the case oi d - wave pairing in the physics of 
copper oxide high - temperature superconductors, more attention will be given to this 
case. 
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Figure 43: Dependence of on the effective width of the pseudogap W in the case of d 

- wave pairing, for several values of impurity scattering rate 'Jo/Tcq: 0-1; 0.18 - 2; 0.64 

- 3. Inverse correlation length /ta=0.2 

In the following we give the results of calculations made for typical values of spectrum 
parameters t' /t = —0.4, fi/t = —1.3, while for correlation length we assume Ka = 0.2. To 
spare space we do not show the results for dimensionless GL - coefficients Ka, Kb, Kq, 
but restrict ourselves to demonstration of more important dependences of all major phys- 
ical characteristics. 

When we are dealing with dependences on the effective width of the pseudogap, all 
characteristics are normalized to their values at T = Tco, while in case of dependences 
on the impurity scattering rate 70 we normalize to the values at T = Tco(W), i.e. at the 
value of the "bare" transition temperature at a given value of W, but in the absence of 
impurity scattering (70 = 0). 

5.5.1 d — wave pairing. 

In Fig. 021 we show the dependence of superconducting transition temperature on the 
effective pseudogap width W for several values of impurity scattering rate. It is seen that 
pseudogap fluctuations lead to signiflcant suppression of superconductivity, and in the 
presence of flnite disorder we always obtain some "critical" value of W, where the value 
of Tc vanishes. This suppression of Tc is naturally due to a partial "dielectrization" of 
electronic spectrum in the vicinity of "hot spots" . H5] . 

Similar dependences for the slope of the upper critical fleld and speciflc heat discon- 
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W/Tco W/Xoo 

Figure 44: Dependence of the slope of the upper critical field and specific heat discon- 
tinuity at the transition on the effective width of the pseudogap W in case of ci - wave 
pairing, for several values of impurity scattering rate 'Jo/Tco- 0-1; 0.18 - 2; 0.64 - 3. 

tinuity at the transition point are shown in Fig. |^ Typically we see fast suppression of 
these characteristics by pseudogap fluctuations. 

Dependence on the value of correlation length of pseudogap fluctuations is more slow 

— in all cases the growth of C, (drop of k) enhances the effect of pseudogap fluctuations. 
We drop appropriate results to spare space. 

In Fig. ^1 we show the dependences of superconducting transition temperature 
on impurity scattering rate 70 for several values of the effective pseudogap width. We 
note that in the presence of pseudogap fluctuations, suppression of Tc by disorder is 
significantly faster, than in their absence {W = 0), when Tc dependence on 70 in case of d 

- wave pairing is described by the standard Abrikosov - Gorkov curve (first obtained for 
s - wave pairing and scattering by magnetic impurities) [TUl lTTj. Similar dependences for 
the slope of Hc2{T) and specific heat discontinuity are shown in Fig. 0^1 Again we see, 
that impurity scattering (disorder) leads to the fast drop of these two characteristics, i.e. 
enhances similar effect of pseudogap fluctuations. 

Similar dependences were obtained in Ref. [72] for the case of pseudogap fluctuations of 
CDW type, when we also obtain recursion procedure for the vertex part wit alternating 
signs. At the same time, in this case we observe some quantitative difference in the 
behavior of all characteristics due to another combinatorics of diagrams. 

Dependences obtained in our model are in qualitative agreement with most of the 
experimental data on superconductivity in the pseudogap region (underdoped region of 
the phase diagram of cuprates). Below we shall show that these results may be used for 
direct modelling of a typical phase diagram of a high - temperature superconductor. 
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Figure 45: Dependence of Tc on impurity scattering rate (disorder) 70 for the case of 0? - 
wave pairing, for several values of the effective pseudogap width: W/Tcq: 0-1; 2.8 - 2; 
5.5 - 3. 




Figure 46: Dependence of the slope of the upper critical field and specific heat discontinu- 
ity at the transition on impurity scattering rate (disorder) 70 in case of d - wave pairing, 
for several values of effective pseudogap width W/Tco'. 0-1; 2.8-2; 5.5-3. 
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Figure 47: Dependence of on effective width of the pseudogap W for the case of s - wave 
pairing, for two values of impurity scattering rate: Jq/Tcq: - 1; 20 - 2. Inverse correlation 
length Ka=0.2. At the insert - characteristic behavior of specific heat discontinuity for 
the same parameters. 

5.5.2 s — wave pairing. 

The case of s - wave pairing is interesting to us mainly with the aim of demonstration of 
major differences from the case of ci - wave case. There are no experimental evidence for 
s - wave superconductivity in systems with pseudogap, though such systems may well be 
discovered in some future. 

Our calculation show that pseudogap fluctuations suppress superconducting transition 
temperature also in this case (Fig. ITTj) . though characteristic scale of these fluctuations, 
necessary for signiflcant suppression of superconductivity here is much larger, than in 
the case of d - wave pairing. This result was obtained in Ref. jT^I, but we must note 
the absence (in the case of Heisenberg (SDW) pseudogap fluctuations) of characteristic 
"plateau" in the dependence of on W, which was obtained for the case of scattering by 
pseudogap fluctuations of CDW type in Ref. j^. On the same scale of W we also observe 
here the suppression of specific heat discontinuity at the transition, which is shown at the 
insert in Fig. |37| 

As to dependence of Tc on impurity scattering rate (disorder), we can see that besides 
relatively small effect of suppression due to [75^ disorder "smearing" of the density of 
states at the Fermi level, even some weak effect of enhancement with the growth of 
7o can be also observed, apparently due to the "filling" the pseudogap in the density of 
states induced by impurity scattering |73j. 

In Fig. HHlwe show the infiuence of impurity scattering (disorder) on the slope of the 
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Figure 48: Dependence of the slope of the upper critical field and specific heat discontinu- 
ity at the transition on impurity scattering rate (disorder) 70 in case of s - wave pairing, 
for two values of pseudogap width: W/Tcq: 0-1; 15-2. 

upper critical field and specific heat discontinuity. Specific heat discontinuity is signifi- 
cantly suppressed by disorder, while the slope of Hc2{T) behaves qualitatively different 
from the case of ci - wave pairing: the growth of disorder leads to the growth of this charac- 
teristic, as in the case of standard theory of "dirty" superconductors |16 , and pseudogap 
fluctuations enhance the slope of Hc2{T). In the absence of pseudogap fluctuations, sim- 
ilar differences between s - wave and d - wave superconductors in behavior of the slope 
of Hc2(T) under disordering were discussed in Ref. [7E] . 

6 Modeling of the phase diagram of cuprates. 

This model for the influence of pseudogap fluctuations on superconductivity allows to 
perform a simple modeling of typical phase diagram of HTSC cuprates [731 E^^^ First 
attempt of such modeling in quite simplified version of our model was undertaken in Ref. 

Basic idea is to identify our parameter W with experimentally observable effective 
width of the pseudogap (temperature of crossover into pseudogap region on the phase 
diagram) Eg ^T*, determined from numerous experiments [211 This characteristic, 
as was already noted above, drops almost linearly with the growth of the concentration 
of doping impurity (concentration of carriers), starting from the values of the order of 
lO^K and vanishing at some critical concentration Xc ~ 0.19. .0.22, slightly higher than 
"optimal" value Xo ~ 0.15. .0.17 P^ I78j. Accordingly, we can assume similar concentration 

-•^^Here we neglect the existence of a narrow region of antiferromagnetism in Mott insulator state, at 
low concentrations of doping impurity, limiting ourselves to a wide region of "bad" metallic state. 
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Figure 49: Model phase diagram for the case of scattering by pseudogap fluctuations of 
CDW - type {d - wave pairing) and "bare" superconducting transition temperature T^q 
linear in carrier concentration. 



dependence of our pseudogap width W{xy^. In this sense we can say that W{x) is 
determined directly from the experiment. Then, the only parameter we need remains 
concentration dependence of the "bare" temperature of superconducting transition Tco(x), 
which would have been existing in the absence of pseudogap fluctuations. Its knowledge 
would allow us to determine the concentration dependence of the observable transition 
temperature Tc{x), solving the equations of our model. Unfortunately, concentration 
dependence Tco{x) is, generally, unknown and is not determined by any known experiment, 
remaining just free fitting parameter of our theory. 

Following Ref. [56j we can assume, that Tco(x) can also be described by linear func- 
tion of X, going to zero at x = 0.3, and choosing the value of Tco{x = 0), which give us 
experimentally observed Tc{x = Xq). Then we can calculate the whole "observed" depen- 
dence Tc{x). Results of such calculation for the case of d - wave pairing and scattering 
by charge (CDW) pseudogap fluctuations [72], using typical W{x) dependence, are shown 
in Fig. 03 We see that even with such arbitrary assumptions the "hot spots" model 
allows to obtain Tc{x) dependence, which is rather close to experimentally observable. 
Similar calculations for Ising like model of spin fluctuations (leading to the absence of 
sign alternation in the recursion procedure for the vertex part [721) show, that reasonable 
values for Tc{x) can be obtained only for unrealistic values of W{x), at least an order of 
magnitude larger than observed in the experiments. 

In the framework of our BCS - like model for the "bare" Tco, an assumption of its 

^^Naturally, such identification can be done up to an unknown proportionality coefficient ^ 1. 
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Figure 50: Model phase diagram for the case of scattering by Heisenberg (SDW) pseudo- 
gap fluctuations {d - wave pairing) and "bare" temperature of superconducting transition 
Tco independent of carrier concentration, with the account of internal disorder linear in 
concentration of doping impurity 7(0;). 

strong concentration dependence seems rather unrealistic^^. Thus, let us assume that the 
value of Tco does not depend on carrier concentration x at all, but take into account the 
fact that introduction of doping impurity inevitably leads to appearance random impurity 
scattering (internal disorder), which can be described by appropriate linear dependence of 
7(0;). Let us assume that it is this growth of disorder, which leads to complete suppression 
of - wave pairing at x = 0.3, in accordance with Abrikosov - Gorkov dependence [7^ 177]. 
Results of our modeling of the phase diagram for the system of the type of La2^xSrxCu04, 
for the case of Heisenberg pseudogap fluctuations, with the account of the abovementioned 
role of internal disorder, are shown in Fig. ^Ul The values of different parameters used in 
this calculation are also shown at the same figure. "Experimental" values of Tc{x), shown 
at this figure (as well as at Fig. by "diamonds" , were obtained from empirical relation 

HHiizni: 

'^"^''^ = 1 _ 82.6(x - Xo)"" (86) 



which gives rather good fit to experimental data on concentration dependence of Tc for a 
number of HTSC cuprates. It is seen, that in the whole underdoped region our model gives 
practically ideal description of "experimental" data with very reasonable values of W{x). 
At the end of overdoped region agreement is less good, but you must take into account. 



"'^''In this approach any dependence of TcO on x can be only due to relatively weak changes of density 
of states at the Fermi level, as chemical potential moves with x. 
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that Eq. (jSHI) also usually does not give very good fit here, and the obvious crudeness 
of our assumptions for underdoped region. Also note, that we have not undertaken any 
special attempts to improve agreement in this region. 

It is interesting to analyze the behavior of superconducting transition temperature Tc 
under additional disordering of our system at different compositions (carrier concentra- 
tions). There are many experimental studies where such an additional disordering 
was introduced by chemical substitutions (impurities) |HI[ IH2] or by irradiation with fast 
neutrons jSD] or electrons [SSI Ell- However, special discussion of the role of this additional 
disordering in the context of the studies of the pseudogap state was made only in Ref. 

In our model, such disordering may be simulated by introduction of additional "im- 
purity" scattering parameter 70, which is just added to the parameter of internal disorder 
7(x). Results of our calculations of superconducting critical temperature for two values 
of this parameter are also shown in Fig. EOl We see, that in complete agreement with 
experiments |H21, introduction of "impurities" (disorder) leads to fast narrowing of the 
superconducting region on the phase diagram. Also in complete accordance with our 
conclusions made above with respect to Fig. E^l and also with experiments [501 E21, su- 
perconductivity suppression by disorder in the underdoped region (pseudogap state) is 
significantly faster, than at optimal composition. We could have expected that introduc- 
tion of "normal" disorder, obviously leading to some suppression of the pseudogap in the 
density of states, could also lead to certain "slowing" down of Tc suppression. However, 
in case of c? - wave pairing, this effect is apparently absent. 

The problem, however, is that in all cases our suppression of by disorder is faster 
than that described by the standard Abrikosov - Gorkov dependence for the case of d 
- wave pairing [77^. At the same time, attempts to fit the results of majority of ex- 
periments on disordering in HTSC cuprates to this dependence show [HH ESI IHll , that 
such suppression actually is significantly slower, than predicted by Abrikosov - Gorkov 
dependence. This, still unsolved, problem remains one of the major problems of the the- 
ory of high - temperature superconductors 16j. One possible solution may be connected 
with consistent description of the role of disorder in superconductors, which belong to the 
crossover region from "large" pairs of BCS theory to "compact" Bosons, appearing in the 
limit of very strong coupling [S3] . Another interesting possibility to explain such "slowing 
down" of Tc suppression is connected with the role of anisotropy of elastic scattering by 
impurities analyzed in Refs. [7^ ISUj . This last effect can be rather easily included in 
our calculational scheme. It is of particular interest in connection with established fact 
of rather strong anisotropy of elastic scattering (with d - wave symmetry), which was 
observed in ARPES experiments on Bi2Sr2CaCu20s+s (SHI- Appropriate scattering rate 
in these experiments changed in the interval of 20 — 60 meV (cf. Fig. El(c)) 3F, which is 
nearly an order of magnitude larger than maximal value of 7(0;), used in our calculations. 
This is an additional evidence of unusual stability of d - wave pairing in cuprates towards 
static disorder. Note that our model for electron self - energy in fact describes similar 
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anisotropy of elastic scattering, corresponding to its growth in the vicinity of "hot spots" , 
yet the effect of "slowing down" of suppression under disordering is not observed. 

Our results show that despite the obvious crudeness of our assumptions, the "hot 
spots" model rather easily leads to a reasonable (sometimes even semiquantitative) de- 
scription of superconductivity region on the phase diagram of HTSC cuprates^^. The 
main shortcoming this approach remains considerable arbitrariness of our "scenario" of 
the origin of concentration dependence of the "bare" superconducting transition temper- 
ature. 

7 Conclusion. 

Our discussion demonstrates the variety of results, which can be obtained in this class of 
models of the pseudogap state. We must stress once again quite simplified nature of our 
approximations, though allowing to obtain what was called "nearly exact" solution. This 
solution is interesting, first of all, from purely theoretical point of view, as rare enough 
example of the problem, when it is possible to sum the complete Feynman perturbation 
theory series (though, in general, certain classes of higher order diagrams are calculated 
approximately). Probably the main theoretical "lesson" to learn is, that the results of 
such complete summation are rather radically different from the results, obtained by one 
or other "partial" summation (e.g. of simple geometrical progression). 

Probably the main shortcoming of our approach is the neglect of dynamics of fluctu- 
ations of short - range order. This is justified, as was shown above, only for high enough 
temperatures, which is bad, e.g. from the point of view of description of superconducting 
state for T < T, 

Another shortcoming, as was mentioned many times [HZl EH]; is our limitation to 
Gaussian fluctuations, which is also justified only for high enough temperatures. In some 
sense it is only technical limitation, though it is quite important for the formal structure 
of our solution. 

We already mentioned that in this class of models we can demonstrate the breaking of 
the fundamental property of self - averaging of superconducting order parameter 
Unfortunately, this can be done only in rather unrealistic case of infinite correlation 
length of pseudogap fluctuations [54j, or in very special model with finite correlation 
length 1^ I59j . and in an oversimplified model with "hot patches" on the fermi surface. 
Qualitative picture emerging from this analysis demonstrates the possibility of formation 
of "superconducting drops" at temperatures higher than mean - field Tc jTHl IMl 1^ - 

^® Above we always assumed that we are dealing with hole doped systems, where concentration depen- 
dence of T*{x) is well established |24l I78j . For electronically doped cuprates similar data are practically 
absent. 

^^Thus, and also to spare space, we have dropped discussion of Gorkov's equations, which can also be 
derived and analyzed |55[l75j . taking into account all higher orders of perturbation theory over pseudogap 
fluctuations. 
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which may have direct relation to the existence of the "low energy" pseudogap (signatures 
of superconducting response at T > T^) and experimentally observed inhomogeneous 
superconductivity in cuprates fri\ I18j. In the "hot spots" model this kind of analysis is 
still to be done. 

At the same time, we have seen that in some cases we can make even direct comparison 
of predictions of these simplified models with the experiments, and the results of such 
comparison are rather promising. In this sense the "hot spots" model can pretend to give 
quasi realistic description of the properties of the pseudogap state in high - temperature 
superconductors. 

The author is grateful to E.Z.Kuchinskii and N.A.Kuleeva (Strigina) for their collab- 
oration during our work on the "hot spots" model. This work was supported in part 
by RFBR grant 02-02-16031 and programs of the Presidium of the Russian Academy 
of Sciences (RAS) "Quantum macrophysics" and of the Division of Physical Sciences of 
the RAS "Strongly correlated electrons in semiconductors, metals, superconductors and 
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